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ABSTRACT 

Context. The AMBRE Project is a collaboration between the European Southern Observatory (ESO) and the Observatoire de la Cote 
d'Azur (OCA) that has been established in order to carry out the determination of stellar atmospheric parameters for the archived 
spectra of four ESO spectrographs. 

Aims. The analysis of the FEROS archived spectra for their stellar parameters (effective temperatures, surface gravities, global metal- 
licities, alpha element to iron ratios and radial velocities) has been completed in the first phase of the AMBRE Project. From the 
complete ESO:FEROS archive dataset that was received, a total of 21551 scientific spectra have been identified, covering the period 
2005 to 2010. These spectra correspond to 6285 stars. 

Methods. The determination of the stellar parameters was carried out using the stellar parameterisation algorithm, MATISSE (MATrix 
Inversion for Spectral SynthEsis), which has been developed at OCA to be used in the analysis of large scale spectroscopic studies 
in galactic archaeology. An analysis pipeline has been constructed that integrates spectral normalisation, cleaning and radial velocity 
correction procedures in order that the FEROS spectra could be analysed automatically with MATISSE to obtain the stellar parame- 
ters. The synthetic grid against which the MATISSE analysis is carried out is currently constrained to parameters of FGKM stars only. 
Results. Stellar atmospheric parameters, effective temperature, surface gravity, metallicity and alpha element abundances, were de- 
termined for 6508 (30.2%) of the FEROS archived spectra (~3087 stars). Radial velocities were determined for 11963 (56%) of the 
archived spectra. 2370 (11%) spectra could not be analysed within the pipeline due to very low signal-to-noise ratios or missing 
spectral orders. 12673 spectra (58.8%) were analysed in the pipeline but their parameters were discarded based on quality criteria 
and error analysis determined within the automated process. The majority of these rejected spectra were found to have broad spectral 
features, as probed both by the direct measurement of the features and cross-correlation function breadths, indicating that they may be 
hot and/or fast rotating stars, which are not considered within the adopted reference synthetic spectra grid. The current configuration 
of the synthetic spectra grid is devoted to slow-rotating FGKM stars. Hence non-standard spectra (binaries, chemically peculiar stars 
etc) that could not be identified may pollute the analysis. 

■ Key words. Methods/data analysis. Astronomical databases/miscellaneous. Stars/fundamental parameters, Techniques/spectroscopic 

' 1. Introduction Spectrometer (RVS) will observe spectra at medium-resolution 

fvj (R^7000 - 1 1,500) which will be used to obtain radial velocities 

^ Astronomy has entered an era of large scale astronomical sur- ^jj targets, as well as to determine stellai- atmospheric pa- 

j, . veys, the scientific goals of which have the potential to consid- rameters and chemical abundances for some ten's of millions of 

erably expand our understanding of the formation and evolution ^^^^ 
of the Universe. In particular current and future large scale spec- 

^ ti-oscopic surveys of the Milky Way will allow asti-onomers to The multi-object instruments at high-resolution (R~20,000; 

: trace in incredible detail the chemical and kinematic history of MIKE on Magellan and FLAMES on the Very Large Telescope), 

- - our Galaxy medium-resolution (R~8,000; AAOmega on the Australian 

These surveys are being undertaken over a range of reso- Astronomical Telescope (AAT)) and low-resolution (R<2,200; 

lutions. Low-resolution surveys (R ^ 2,000 to 8,000) are now ^'"^'^ ^ulti Object Spectrograph (FMOS) on the Subaru 

widespre ad, for exam ple RAVE (ISteinmetz et all IM and Telescope) are also being used to address science goals specific 

SEGUE (lYannvetal.i r2009). while high-resolu tion surveys are '° ^^^^ of galactic archaeology One key instrument that is 

a more recent endeavour such as APOGEE (-MaiewskietalJ ^^'"^ ^^'^^ ^"'^ '^^'^^"^'^ pnmanly for galac- 

l2007h ) and the Gaia-ESO Survey (Pl.s: Gerry Gilmore & Sofia archaeology research is the H^fh Effi ciency and R egglution 

Randich). The science goals of these ground-based surveys will Multi-Element Spectrograph (HERMES | Barden et al. || 2010D on 

significantly contribute to the studies of galactic archaeology ^ 

as well as provide complementary information to the upcom- This ongoing accumulation of large spectral datasets from 

ing astrometric survey, the European Space Agency (ESA) Gaia surveys and individual spectrographs has compelled the devel- 

Mission. The ESA Gaia satellite will observe approximately opment of automated stellar parameterisation algorithms that 

a billion stars in the Milky Way for which distances will be can reliably and effectively analyse every spectrum. The stellar 

determined to milliarcsecond accuracies. Its Radial Velocity parameterisation algorithm, MATISSE, has been developed at 
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the O bservatoire de la Cote d'Azur (OCA) (iRecio-Blanco et a"D 
I2006 h to be included in the automated pipeline for the analysis 
and parameterisation of the Gaia-RVS stellar spectra. MATISSE 
has also been developed as a standalone Java app lication for use 
in a wide variety of projects (See for instance, 'G azzano et alj 
120 lot iKordopatis eTal] l2011h . The AMBRE Project team at 
OCA, which oversees the development of MATISSE, is con- 
nected to the Gaia Data Processing Consortium (DPAC) under 
the Generalized Stellar Parametrizer-spectroscopy (GSP-spec) 
Top Level Work Package which is overseen by Coordination 
Unit 8 (CU8). 

With such large datasets soon to be available a crucial as- 
pect of the analysis is to be able to compare the results from 
the different surveys. To do, this a comprehensive set of stan- 
dard objects is required that can be used to calibrate the datasets. 
Preparation for this is in the form of the development of spec- 
tral libraries. These are datasets of spectra with homogeneously 
determined characteristics, such as stellar parameters, which can 
be used as calibration stars for these surveys. In particular stan- 
dard star lists are being developed for Gaia for both the radial 
velo city (Crifo et al. 2010) and stellar parameter measurements 
JSoubiran et al.ll2010l) . 

The work carried out for the AMBRE Project in effect con- 
verts the spectra in the European Southern Observatory (ESO) 
archive into a comprehensive spectral library of homogeneously 
determined stellar parameters. There are three primary objec- 
tives for the AMBRE Project: 

1. To rigorously test MATISSE on large spectral datasets over 
a range of wavelengths and resolutions, including those for 
the Gaia RVS, 

2. To provide ESO with a database of stellar temperatures, 
gravities, metallicities, alpha to iron ratios and radial veloci- 
ties for the associated archived spectra that will then be made 
available to the international scientific community via the 
ESO Archive, 

3. To create a chemical map of the Galaxy from the combined 
ESO archived sample upon which stellar and galactic for- 
mation and evolution archaeological analysis can be carried 
out. 

The first phase of the AMBRE Project was the analysis of 
the FEROS archived spectra. From the complete ESO:FEROS 
archive dataset that was delivered to OCA, a total of 21551 sci- 
entific spectra have been identified, covering the period 2005 to 
2010. These spectra correspond to 6285 different stars based on 
a coordinate matching calculation with a radius of 10". 

The structure of this paper is as follows: Section|2]introduces 
the AMBRE Project; Section [3] introduces the MATISSE algo- 
rithm and the synthetic spectra grid; Section|4]describes the anal- 
ysis pipeline that has been built around MATISSE for the analy- 
sis of the archived spectra; Section|5]discusses the internal errors 
that have been calculated for the pipeline; Section |6] discusses 
the external errors analysis based on a reference sample of stars; 
Section|2]presents the application of key rejection criteria to the 
spectral dataset; Section [8]presents the stellar parameter results 
for FEROS and finally Section |9]concludes the paper 

2. The AMBRE Project 

Under a contract with ESO the archived spectra of four 
ESO spectrographs are being analysed using the auto- 
mated parameterisation programme MATISSE in a project 
overseen by the AMBRE Project team. Table [1] lists the 



Table 1. Details of the four ESO spectrographs and their publicly 
available archived spectra sample that are part of the AMBRE 
Project. 



Spectrograph 


Resolution 


A (nm) 


No. spectra 


FEROS 


48,000 


350- 


-920 


21,551 


HARPS 


115,000 


378- 


-690 


126,694 


UVES 


20,000-110,000 


300- 


1100 


78,593 


FLAMES/GIRAFFE 


5,600-46,000 


370- 


-900 


> 100,000 



main characteristics of the four spectrographs: the Fiberfed 
Extended Range Optical Spectrograph (FEROS) dKaufer et al.l 
Il999h : th e High Accuracy Radial velocity Planet Searcher 
(HARPS) ('Mayor et al."2003'); the Ultraviolet and Visual Echelle 
Spectrograph (Dekkeretal. 2000); and the Fibre Large Array 
Multi Elemen t Spectrograph/GIRAFFE (FLAMES/GIRAFFE) 
dPasquini et al..i2002) . The atmospheric stellar parameters of ef- 
fective temperature (Tes), surface gravity (logg), global metal- 
licity ([M/H]), a element to iron ratio ([c/Fe]) and radial veloc- 
ity (Vrad) will be derived for each of the archived stellar spectra. 
These will be delivered to ESO for inclusion in the ESO database 
and then made available to the astronomical community via the 
ESO archive. It is intended that the availability of stellar param- 
eters for each spectra will encourage further use of the archived 
spectra. Previously unconsidered samples can be found through 
searches on the parameters, for examples extracting all the spec- 
tra in a particular metallicity range, or with very similar temper- 
atures and gravities. 

The analysis of the archived spectra of these four spectro- 
graphs presents a unique opportunity to test the performance 
of MATISSE on large datasets of real spectra. In particular 
key instrument configurations of FLAMES/GIRAFFE cover the 
Gaia RVS wavelength domain and resolutions. Rigorous test- 
ing of MATISSE is necessary in order to optimise its perfor- 
mance in the Gaia-RVS analysis pipeline that is being compiled 
at the Centre National d'Etudes Spatiales (CNES). As such the 
AMBRE Project has been formally designated as a sub-work 
package under GSP-spec. The stars analysed by AMBRE will 
also be available for use as standard or calibration stars for 
the Gaia-ESO survey, and as secondary standards for the Gaia 
Mission. 

The first phase of the AMBRE Project, the analysis of 
the FEROS archived spectra, is now complete. FEROS is 
a state-of-the-art bench-mounted high resolution spectrograph 
that was built by a consortiurrQ of four astronomical insti- 
tutes: Landessternwarte Heidelberg, Astronomical Observatory 
Copenhagen, Institut d' Astrophysique de Paris and Observatoire 
de Paris/Meudon. In 1998 FEROS was installed and commis- 
sioned on the ESO 1.52 m telescope at La Silla, Chile. In 2002 
FEROS was moved to the MPG/ESO 2.2 m telescope where it 
permanently residefl The high resolution and full coverage of 
the optical domain means that FEROS is suitable for a wide 
range of astronomical projects from radial velocities and vari- 
ability studies to chemical abundances in stellar populations. 

In 2009-10 the FEROS spectra from October 2005 to 
December 2009 were reanalysed by the ESO archivists using 
an improved reduction pipeline. These reduced spectra were the 
sample that was delivered to OCA for analysis in AMBRE. 

Figure [1] shows a histogram of the number of FEROS obser- 
vations per year that were received by OCA, the number for 2005 

' http://www.lsw.uni-heidelberg.de/projects/instrumentation/Feros/ 
" http://www.eso.org/sci/facilities/lasilla/instruments/feros/index.html 
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Fig. 1. Number of FEROS observations per year for anal- 
ysis in AMBRE. The 2005 observations begin in October 
The number of distinct objects observed and the number 
with repeated observations are as indicated in the key. 



representing only 3 months of observations. The number counts 
for the number of distinct objects that have been observed by 
FEROS in this timeframe and the number with repeated obser- 
vations are also shown. 

The analysis procedure and results for these FEROS spectra 
are presented in this paper This analysis has been the testbed for 
producing many of the tools that will also be used in the analysis 
of the UVES, HARPS and FLAMES/GIRAFFE archived spec- 
tra. These tools have been integrated into an analysis pipeline 
that feeds the processed spectra into MATISSE for derivation of 
the spectral parameters. 

3. MATISSE & the synthetic spectra grid 

MATISSE (MATrix Inversion for Spectral SynthEsis) is an auto- 
mated stellar parameterisation algorithm based on a local multi- 
linear regression method. It derives stellar parameters (9 - T^ff, 
log g, [M/H], individual chemical abundances) by the projection 
of an input observed spectrum on a vector function Bg{A), 



9i = J] BeiA) ■ OiiA) 



(1) 



with 9 being the derived value. The Be(A) vector function is an 
optimal linear combination of theoretical spectra, S i{A), calcu- 
lated from a synthetic spectra grid. 



Be,(A)^Y, aijSj(A) 



(2) 



Key features in the observed spectrum due to a particular 
9 are reflected in the corresponding Bg{A) vector ind icating the 
particular regions which are sensitive to 9 dRecio-Blanco et alj 
l2006HBiiaoui et al.ll2008h 1. 

MATISSE also generates synthetic spectra interpolated to 
each set of output stellar parameters. For each set of parameters, 
the synthetic spectrum on the grid, SoiA), with the closest stellar 
parameters is identified along with the associated Be{A) function. 
The Bg{A) functions located locally about this point are used to 



estimate the variations in the flux between the synthetic spectra 
at the corresponding grid points and the final interpolated syn- 
theti c spectrum correspond ing to the required stellar parameters 
(9ic) (iKordopatis et al 

This interpolated spectrum is then used to calculate a be- 
tween the interpolated and input normalised spectrum. This pro- 
vides a measure of the goodness of the fit of the derived stellar 
parameters to the observed spectrum. Stellar parameters, with 
logx^, and corresponding interpolated synthetic spectra are gen- 
erated for every spectrum that is analysed in MATISSE. 

3.1. Grid of synthetic spectra for tlie classification algoritlim 

Within the AMBRE Project, the adopted procedure for the auto- 
matic classification of stellar spectra relies on a library of refer- 
ence spectra in order to derive their atmospheric parameters and 
their chemical abundances. Due to the lack of a library of ob- 
served spectra that covers a large range of atmospheric parame- 
ters and chemical abundances over a very large spectral domain 
and resolution, the only solution was to compute large grids of 
synthetic spectra. 

For the AMBRE application, the computed grid of theo- 
retical stellar spectra has to cover the whole optical spectral 
range at very high resolution in order to be used for the anal- 
ysis of the majority of the spectroscopic data. For that reason, 
we have computed a synthetic spectra grid covering the wave- 
length range between 300 and 1 200 n m with a wavelength step 
of 0.001 nm (900,000 pixels in total) jde Lavernv et al.ll2012l) . 
Since this project is mostly devoted to the analysis of FGKM 
stars, this grid is based on the latest generati on of MARCS model 
atmospheres presented in lGustafsson et all (.2008 ). An extension 
to hotte r effective tem peratures with Kurucz stellar atmosphere 
models (iKurucdl979h is planned for the near future. The consid- 
ered parameters of the spectral grid are the effective temperature 
(Teff in K), the stellar surface gravity (logg in dex), the mean 
metallicity ([M/H] in dex) and the enrichment in a-elements 
with respect to iron ([a/Fe] in dex). 

The metallicity definition here is different to the classical 
use of [Fe/H] to denote the metallicity of a star. For [Fe/H] the 
metallicity is defined using the derived abundances from Fe lines 
only. MATISSE provides the opportunity to use all the available 
metal lines (those atoms heavier than He) to define a metallicity 
providing a global metallicity designated as [M/H]. The oppor- 
tunity also exists with MATISSE to derive a global a element 
abundance using as many a element spectral lines as possible, 
where we assumed that in the generation of the synthetic spec- 
tra the abundances of the different a elements vary in lockstep. 
Throughout the AMBRE Project, the following chemical species 
are assumed to be a-elements: O, Ne, Mg, Si, S, Ar, Ca and Ti, 
although for any selected wavelength region spectral features for 
all of these elements may not necessarily be present. 

From the selected MARCS model atmospheres, the syn- 
thetic spectra were c omputed with the turbospectrum code 

varez & Pled 1 19981 and further improvements by Plez) in 
plane-parallel and spherical geometry assuming hydrostatic and 
local thermodynamic equilibrium. Atomic lines have been re- 
covered from the Vienna Atomic Line Database (in August 
2009: iKupka et alJI 1999 ). The molecular line list has been pro- 
vided by B. Plez. It includes transitions from ZrO, TiO, VO, OH, 
CN, C2, CH, SiH, CaH MgH and FeH with the ir corresponding 
isotopic variations (see lGustafsson et al.l (l2008l) . for a list of ref- 
erences). 

This grid covers the following ranges of atmospheric pa- 
rameters: Tgff between 2 500 K and 8 000 K, \ogg from -0.5 
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to +5.5 dex, and [Fe/H] from -5.0 to +1.0 dex, although 
not all combinations of the parameters are available within 
the grid. The selected MARCS models have [afFe]=Q.Q for 
[M/H] > 0.0, [a/Fe]=+0.4 for [M/H] < -1.0 and, in between, 
[a/Fe]=-0.25x[M/H]. For the spectra computation from each 
of these MARCS models, we considered an [a/Fe] enrichment 
from -0.4 to +0.4 dex with respect to the canonical values that 
correspond to the original abundances of the MARCS models. 
The final AMBRE synth etic spectra grid consis ts of 16783 flux 
normalized spectra (See Ide Lavernv et al.ll2012[ for a complete 
description of this grid). 

The microturbulence (^) is not a free parameter in this syn- 
thetic spectra grid. For the atmospheric models with high log g 
(+3.5 < \ogg < +5.5) ^ is set at 1.0 kms"'. For low \ogg 
(logg < +3.0) ^ is set at 2.0 kms"', being typical values for 
dwarfs and giants respectively. These values reflect the MARCS 
model atmospheres config urations, and further de tails on the 
model selection is given in Ide Lavernv et all (|2012|) . Due to the 
microturbulence being hardwired into the synthetic grid it was 
not possible to carry out tests on the effects of variations in 
^ on the stellar parameter determination across the whole syn- 
thetic spectra grid. However the effects of changes in ^ are most 
prominent for strong lines. By using the global metallicity [M/H] 
rather than [Fe/H] we expect the contribution from strong lines 
on the [M/H] to be negligible due to the significantly larger quan- 
tity of non-^ sensitive small metallic lines. 

Section|6]describes the extensive testing of the pipeline that 
was carried out using observed spectra from which an exter- 
nal error on the resulting parameters is derived by compari- 
son to literature values. To gauge the effect of microturbulence 
on the derived stellar para meters a key sample of 66 Main 
Sequence (MS) stars (Ben sbv et al.l l2003l) . and two key sam- 
ples of 16 Red Gian t Branch (RGB) stars (lMcWilliamlll990t 
iHekker & Melendezlllool were considered from within the test 
sample. 

The most noticeable effect within the MS sample was an 
underestimation of the Teff for the hotter MS stars, for which 
the constant ^ of 1.0 kms"' underestimates the accepted value. 
Hence at Teff ~ 6500 K, where ~ +0.8kms"' the effect was 
Areff ~ -100 K. This is well within the derived T^ff external 
error of 120 K (see Section|6]l. 

The most noticeably effect in the RGB sample was on the 
derivation of logg. At the base of the RGB, where the constant 
^ of 2.0 kms"' is an overestimation, the derived logg was ob- 
served to be overestimated by ~ 0.3 dex for ~ +0.7kms"'. 
At the RGB tip, where ^ = 2.0 kms ' is an underestima- 
tion, logg was observed to be underestimated by ~ 0.3 dex for 

A^ 0.7 kms"'. These variations are within the derived logg 

external error of 0.37 dex. Hence based on these samples, vari- 
ations of ^ from the assumed constant values do have an effect 
on the derived parameters that varies in magnitude with stellar 
evolutionary stage. However the external error derived from the 
entire reference sample takes account of them in a global sense 
(see Section|6]i. 

3.2. AMBRE:FEROS subghd 

The FEROS spectra cover most of the optical wavelengths 
but only a subset of these wavelengths was required for the 
MATISSE analysis. In order to caiiy out the MATISSE train- 
ing phase, where the Be{A) vectors are generated, for the FEROS 
analysis the optimum resolution, wavelength regions, and sam- 
pling of the FEROS spectra were determined. The full AMBRE 
synthetic spectra grid was then adapted to the same specifica- 



tions creating the AMBRE:FEROS subgrid that was used in the 
training phase to generate the AMBRE:FEROS Bg(A) vectors. 

In this analysis, only Bg{A) functions computed from 
the direct numerical in version of the correlation matrix (see 
iRecio-Blanco et 31112006) were considered. These functions are 
thus not optimized for very low SNR spectra. However, due to 
the large number of spectral features available in the selected 
spectral domains, it was not necessary to use the approximated 
Bg(A ) functions that are comp uted with the Landweber algorithm 
as in lKordopatis et al 1 (1201 ij). As will be shown in Section|5]the 
estimated internal errors based on the AMBRE:FEROS grid are 
indeed already very small (see Figure [TOt. 

3.2.1. Selection of wavelength regions 

A detailed analysis of the full wavelength range of the FEROS 
spectra was undertaken in order to select wavelength regions that 
would provide the greatest amount of information for each of the 
derived parameters while also minimising the number of pixels 
in order to reduce computing time. 

FEROS disperses the spectra into 39 orders of varying wave- 
length interval from ~3500 A to ~9200 A at high resolution 
(R~48,000). The optimum spectral regions for the determination 
of the stellar parameters were selected taking into account the 
need to avoid excessive computing time and to avoid low spec- 
tral information region. To this end, the signal-to-noise (SNR) 
per pixel (0.03 or 0.06 A per pixel for FEROS) as a function of 
the wavelength was derived for each archived spectra by deter- 
mining the SNR profile of each spectral order and thus providing 
a SNR profile over the entire spectral domain. Wavelengths re- 
gions with low SNR were rejected, typically the start and end of 
each order The wavelengths regions affected by sky absorption 
and telluric features were also rejected, and also regions where 
continuum placement proved too difficult owing to wide spec- 
tral features in the region. Some key spectral features, such as 
H„ and the Ca 11 H & K lines, were also discarded as these fea- 
tures were found to be poorly synthesised based on our current 
understanding of stellar atmospheres as well as being difficult to 
normalise automatically. 

The remaining spectral regions were then investigated for 
their intrinsic sensitivity to three of the four stellar parameters 
to be determined by MATISSE: T^ff, logg; and [M/H]. This was 
carried out using preliminary Bg{A) vectors that covered the full 
optical domain at very high resolution generated using a reduced 
synthetic grid. Due to the underlying equation for calculating 
the Bg{A) vectors, the high resolution Bg(A) vectors could not be 
simply degraded to lower resolutions to investigate the sensitiv- 
ity, as this would not accurately reflect the information at that 
resolution. However the high resolution Bg{A) vectors are useful 
to give a general sense of the sensitivity. 

For stellar parameters at: 

1. reff=6000 K, logg=4.0 dex, [M/H]=0.0 dex (Sun); 

2. reft=6000 K, logg=2.0 dex, [M/H] =0.0 dex; 

3. refF=4500 K, logg=2.0 dex, [M/H]=-0.5 dex (Arcturus); 

4. reff=4500 K, logg=4.0 dex, [M/H] =-0.5 dex, 

the Bg(A) vectors at the same T^ff, logg, and [M/H] have been 
analysed for their extrema distribution. This selection explores 
the sensitivity of both dwarfs and giants that are similar to the 
two standard stars, the Sun and Arcturus. The panels in Figure|2] 
show how the spectral domain changes in sensitivity with the 
different stellar parameters. The standard deviation (cr) for each 
Bg{A) vector was calculated as the spread of the pixel values in 
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Fig. 2. Bg(A) Sensitivity: Number of extrema above 2 cr in 
~200 A bins against wavelength for four sets of stellar parame- 
ters: the Sun; T^s = 6000 K, logg = 2, [M/H] = 0.0; Arctums; 
and Teff = 4500 K, logg = 4, [M/H] = -0.5. The Be{A) vec- 
tors from left to right are: 6 - T^ff, logg and [M/H]. The grey 
stripes represent the final wavelength regions selected for the 
AMBRE:FEROS analysis (see Table©. 
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Fig. 3. As for Figure|2]but for the corresponding BgiA) vectors of 
the AMBRE:FEROS synthetic spectra grid. The y-axis is scaled 
down compared to Figure |2] 



the Bg(A) vector Then for each Bg(A) vector the pixels outside 
2 (T were binned in ~200 A bins over the wavelength range from 
~4000 A to ~9500 A. Hence the bars in Figure |2]show the num- 
ber of pixels per bin which have a high sensitivity to the respec- 
tive 0. Clearly the bluer wavelengths show the greatest sensitivity 
to all three parameters, reflecting the greater quantity of spectral 
features in the blue. 

Table 2. FEROS echelle order, starting wavelength and finishing 
wavelength for each region used in the analysis of the FEROS 
spectra. 



Region A min A max Region A min A max 





A 


A 




A 


A 


1 


6660 


6730 


10 


4585 


4670 


2 


5655 


5725 


11 


4505 


4580 


3 


5390 


5430 


12 


4405 


4492 


4 


5260 


5320 


13 


4300 


4405 


5 


5120 


5190 


14 


4245 


4300 


6 


4990 


5055 


15 


4140 


4220 


7 


4910 


4960 


16 


4075 


4140 


8 


4800 


4890 


17 


4019 


4075 


9 


4674 


4740 









This sensitivity to d for different ranges in was used to se- 
lect the wavelengths regions in the FEROS spectra to be used 
in the AMBRE analysis. As the greatest sensitivity was located 
towards the blue, regions were selected, where possible, to cap- 
ture this sensitivity. The grey regions in Figure |2]represent the 
final wavelength regions selected for AMBRE:FEROS which 
are listed in Table |2] The wavelengths include key spectral fea- 
tures such as the magnesium triplet at ~5160 A, H^ (4861 A), Hy 
(4340 A), H^ (4101 A), a CH bandhead (4305 A) and a CN band- 
head (4142 A) as well as the spectral features for many atomic 
elements. This wealth of information provided sufficient sensi- 



tivity over the range of required stellar parameters in anticipa- 
tion of the range of spectral types that were potentially included 
within the FEROS archived sample. 

For confirmation of the Be(A) sensitivity. Figure |3] was gen- 
erated by the same process as Figure |2]but using the correspond- 
ing Bg{A) vectors for the final AMBRE:FEROS grid. Therefore 
Figure |3]reflects the sensitivity of the Bg(A) vectors at the actual 
resolution of the AMBRE:FEROS analysis. Due to the reduced 
number of points in sampling and resolution the 2 cr limit is 
significantly reduced. While some regions show there are subtle 
differences in the sensitivity of some regions, overall the distri- 
bution of greater sensitivity in the blue is in good agreement with 
the sensitivity at high resolution. 

3.2.2. Convolution & sampling of synthetic spectra grid 

At the high resolution of R~48,000, the FEROS spectra have 
typical wavelength sampling (AA) of 0.03 A (1x1 binning) or 
0.06 A (1x2 binning). The final selected wavelengths corre- 
sponded to a total of ~ 1500 A and at this sampling this translated 
to a (maximum) total of ~50,000 pixels. Creating the FEROS 
synthetic spectra subgrid of ~ 16,000 spectra at this resolution 
and sampling would result in excessively large memory and 
computing requirements when creating the Bg(A) vectors in the 
training phase, and also when loading these vectors during the 
MATISSSE analysis. It was necessary to reduce the resolution 
and sampling to a more reasonable pixel total, but without sacri- 
ficing the key spectral information. The most reasonable number 
of pixels per spectra was determined to be ~15,000 in order to 
ease the use of the available computing power 

By optimising the memory and spectral information require- 
ments, the final specifications for the FEROS synthetic spectra 
subgrid was determined to be 11890 pixels with R~15,000 at 
A ~ 4500 A with a sampling of 0. 1 A over ~1500 A. Hence 
the full AMBRE synthetic spectra grid (Ide Lavemv et al.ll2012l) 
was then sliced to the wavelengths specified in Table |2l con- 
volved by a Gaussian kernel with constant cr and resampled to 
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the optimised FEROS pixels. The training phase was then car- 
ried out whereby the Bg{A) vectors were generated for use in the 
MATISSE analysis of the FEROS archived spectra. 

4. FEROS analysis pipeline 

A complex analysis pipeline (written in shell, Java and IDL) 
has been built to wavelength slice, radial velocity correct, nor- 
malise and convolve the FEROS spectra and then feed them 
into MATISSE for the determination of their stellar parameters. 
Figure |4] shows a flowchart of the key stages in the analysis 
pipeline. The pipeline has been developed so that it can be eas- 
ily adapted to the remaining three instrument datasets of ESO 
archived spectra. The following sections outline the procedures 
carried out for each of the key stages in the pipeline. 



FEROS Analysis Pipeline 



■* V^,^, FWHM Of CCF 

Error estimate 
■* Preferred -mask parametei 



Wavelength selection 
Basic Normalisation 
Cosmic ray removal 
Determine Signal-to-Noise 
Assumed FWHM 




Modified Spectral Processing B 
Convolution with measured FWHM 
Normalisation using synthetic spectra 
■* Re-Determine Signal-to-Noise 



-> Stellar parameters 
-» Tests 



Final Stellar Parameters 

Final Synthetic Spectra 
Final Radial Velocity 



Fig. 4. The FEROS Analysis Pipeline: The key stages are dis- 
played in order of analysis. Spectral Processing A and Spectral 
Processing B carry out testing of spectra quality and preliminary 
parameter determination, including calculation of the radial ve- 
locity and spectral FWHM. In Spectral Processing C robust it- 
erative procedures are carried out resulting in the final stellar 
parameters and normalised spectra. See text for detailed expla- 
nation of each stage. 



4.1. Spectral Processing A: Spectra processing for radial 
velocity determination 

Spectral Processing A (SPA) is the initial stage that prepares 
the FEROS spectra in order to determine the radial velocities. 
The procedures used here are also used in the later stages of the 
pipeline. Key flags were defined within the pipeline to be at- 
tached to any spectra that satisfied the conditions of the flag. 
The flags were for 'Faulty Spectra' (e.g. missing spectral or- 
ders), 'Extreme Emission Features' (spectral emission features 
in the blue with width greater than 50 A), 'Excessive Noise' 
(extremely noisy spectra with the number of negative flux pix- 
els > 25%) and 'Poor Normalisation' (spectra with a total flux 
above 1 .2 > 10%). A spectrum may satisfy multiple flags, partic- 
ularly in the case of excessive noise which naturally led to poor 
normalisation. Such types of spectra were unlikely to produce 
valid results in the MATISSE analysis and as such were rejected 
from the full analysis. Flags were also implemented to iden- 
tify 'Large Emission Features' (width greater than 1.2 A) and 



the presence of odd features that were potentially 'Instrumental 
Relics'. However these flags alone were not enough to reject a 
spectrum from the full analysis. All of the quality flags were ac- 
counted for in the generation of the final dataset of results for 
ESO, which is explained in Section[8] 

The first reduction process extracts the required wavelengths 
(Table |2| from the spectra and then normalises them. However 
no convolution is carried out so the original resolution is main- 
tained. A hot star spectrum (almost featureless) in the FEROS 
sample was identified and used to determine a continuum pro- 
file for each wavelength region. This profile is divided out of the 
spectra removing the residual instrument profile and performing 
a first rough normalisation to place the continuum near unity. At 
this stage the spectra must necessarily be considered as unknown 
in stellar parameters and hence the normalisation process must 
be canied out assuming that key pre-selected continuum regions 
need to be normalised to unity in intensity. This assumption may 
not in fact correctly normalise each spectrum. However in these 
initial stages it is the best, albeit crude, estimate available. 

Up to three continuum regions for each wavelength region 
were pre-selected by an investigation of synthetic spectra and a 
sample of FEROS spectra covering a range of stellar parame- 
ters. For the normalisation, a simple linear profile is constructed 
for each wavelength region by determination of the mean signal 
in the respective continuum regions. This profile is divided out 
of the wavelength region and then the region undergoes clean- 
ing for cosmic rays. The region is scanned for emission features 
greater than a factor of 1.5 above the noise on the continuum 
(unity). These features are tested for width and those less than 
1 .2 A (40 pixels) are cleaned to unity. Features larger than 1 .2 A 
(40 pixels) in width are flagged as large emission features and are 
not cleaned. This limit of 1 .2 A (40 pixels) was determined by 
an examination of many examples of cosmic rays in the FEROS 
spectra to get a sense of the number of pixels over which they 
could be dispersed. 

Each wavelength region is treated individually for each spec- 
trum, and the regions are then combined into one vector This 
merged spectrum is tested for the goodness of the normalisation 
and flagged if the normalisation is poor The resulting spectrum 
is used in the determination of the radial velocity. At this stage 
the original resolution is retained. None of the FEROS spectra 
were rejected at this stage in case later testing and processing re- 
covered a previously rejected spectra (i.e. mis-calculated radial 
velocity is corrected). 



4.2. Radial velocity determination 

A completely automatic radial velocity {V rad) programme has 
been established that can analyse spectra across a very wide 
range of stellar parameters and is based on a cross-correlation 
algorithm that compares the normalised obser ved spectrum to bi - 
nary masks (C. Melo, private communication. iMelo et al.ll200lh . 
To complement this programme a procedure was developed at 
OCA that computes a set of binary masks to match the obser- 
vations in terms of wavelength range and resolution, and are 
used as input to the radial velocity programme. Masks covering 
the parameter space of the AMBRE:FEROS analysis were com- 
puted using synthetic spectra taken from the high resolution, full 
optical domain, AMBRE synthetic spectra grid (see SectionjTTJ. 
Six masks for stars hotter that 8000 K were also compu ted us- 
ing synthetic spectra from the POLLUX database (Palaci os et al.l 
2010). For the AMBRE:FEROS analysis 51 masks were com- 
puted at the stellar parameters given in Table |3] Also listed are 5 
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standard masks for dwarf stars that were supplied with the radial 
velocity programme (C. Melo, private communication) and used 
in the AMBRE:FEROS analysis. 

For each spectrum a cross-correlation function (CCF) rang- 
ing from -500 to -1-500 kms"' was calculated for each of the 56 
masks. The CCF was calculated in radial velocity steps (AYi-ad) 
based on the specified wavelength step (AA) of the spectrum. 
For spectra with AA = 0.03A the CCF step was calculated to 
be AW i-ad ~ 1-8 kms"', and for AA = 0.06A the CCF step was 
^rad ~ 3.6 kms"'. For each CCF two separate gaussian fits 
were made in order to determine the minimum of the profile, 
and hence the radial velocity. The first fit used the majority of the 
profile, while the second fit used a section of the profile centred 
near the mininum with a width equal to the full-width-at-half- 
maximum (FWHM) of the profile. 




= -46.9±0.1kms' 



Fig. 5. CCFs produced for calculating radial velocities, a) Noisy 
CCF with an asymmetric profile, b) Core of profile showing 
skewed full profile gaussian fit (black). The secondary gaussian 
fit (red) provides a better estimate, c) CCF of the spectroscopic 
binary, HD 135728, showing two profiles. 



Figure |5^. shows an example of a noisy CCF although the 
key profile is distinct. However, as seen in Figure ISji., the fit to 
the full profile (black) provides a minimum skewed away from 
the actual minimum due to asymmetries of the full profile. The 
use of the central core of the profile provides a more accurate fit 
(red). Figure |5];. shows an example of a CCF for which two pro- 
files are observed. This is the spectroscopic binary HD135728. 
The pipeline selects the most prominent of the two profiles in the 
determination of the radial velocity. The CCF seems the most 
likely tool with which to be able to identify the spectra of spec- 
troscopic binaries. However at this stage we have not developed 
such a routine and so binaries are not specifically detected in the 
pipeline. The stellar parameters for such spectra are likely to be 
poorly determined and so will be at least identified as having 
large associated errors from the radial velocity and stellar pa- 
rameter analysis. Quality control flags, such as the x^i may also 
be used to help identify spectroscopic binaries. 

Figure |6] shows the CCF corresponding to the best fit mask 
for three FEROS spectra at different stellar parameters. The full 
CCF for each spectrum is shown in the first row, while the sec- 
ond row shows more clearly the profile used to calculate the ra- 
dial velocity. The full CCF can show varying degrees of noise 
and gradient depending on the quality of the spectra as is shown 
in each of these three examples. The respective profiles used 
to determine the radial velocity reflect the characteristics of the 
different spectra. The two cooler stars have very prominent in- 
verted peaks, reflecting the numerous spectral lines available 
with which to make the cross-correlation. The hotter star, with 
fewer spectral lines, has a much less pronounced profile. These 
differences can be quantified by the measurement of the contrast 
of the profile. 



The contrast was defined here as the amplitude of the pro- 
file multiplied by 100 (hence a percentage) where the amplitude 
of the profile is the length of the vertical red line as shown in 
each example profile in Figure |6] The amplitude of the CCF, 
the continuum placement of the CCF and their associated errors 
are calculated automatically when using the GAUSSFIT routine 
in IDL. The sign of the amplitude does change depending on 
whether the CCF has an absorption or an emission profile. The 
expected profile in this analysis is an absorption profile. 

The two conditions that were required for a profile to have 
been well fitted by a gausssian were first: the contrast must be 
less than zero ensuring an absorption profile; and second, the 
error on the amplitude must be less than 20% of the ampli- 
tude ensuring that the CCF profile dominates above the noise. 
CCFs that satisfied these conditions were considered to be well- 
defined. Finally the error on the radial velocity (cryi-ad) was cal- 
culated using the prescription outlined in lTonrv & Davisl(ll97^ . 
This prescription makes use of the relative heights of the primary 
and secondary peaks in the CCF profile hence crvrad reflects the 
goodness of the definition of the primary peak. 

The radial velocity programme carried out the calculations 
for each of the 56 masks, resulting in 56 determinations of the 
radial velocity for each spectrum. The errors derived from the 
quality of the gaussian fit and crvrad were used to select the best 
fit radial velocity. Typically the majority of these determinations 
were in good agreement at the correct radial velocity. However 
for masks with parameters far from the true parameters of the 
star, a radial velocity would be determined that was incorrect 
but sometimes had sufficiently small errors such that it would 
be incorrectly selected as the best fit, if the smallest error was 
the only selection criteria. To avoid this a binning procedure was 
implemented as follows that discarded such outliers. 

For a single spectrum, all of the Wrad determinations with 
well-defined CCFs were binned by N rod into a maximum of 3 
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Fig. 6. Examples of cross-correlation functions of the best fit 
masks for three of the FEROS spectra. The first row shows the 
full CCFs while the second row shows the profiles from which 
the radial velocities were determined. The final AMBRE:FEROS 
stellar parameters (Teff, log g, [M/H], [a/Fe]) and radial veloci- 
ties are stated for each spectrum. The red lines are used to calcu- 
late the contrast and FWHM for each profile. 
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Table 3. List of the stellar parameters for the binary masks used in the AMBRE:FEROS radial velocity determination that were 
built from the full AMBRE synthetic grid (45 masks) and the POLLUX database (6 masks). Five standard masks for dwarf stars 
provided by C. Melo are also listed. 
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bins in correlation with the overall spread of the V rad values. The 
bin size was set in each instance by the difference in the maxi- 
mum and minimum V rad values divided by the number of bins. 
The bin with the largest number of members was retained and 
the outlying bins were discarded. The N rad within accepted bin 
with the smallest crvrad was selected as the final N rad- In the case 
of an equal number of members within 2 or more bins, the bin 
with the smallest spread in V rad values (so best agreement) was 
selected. Typically, most of the radial velocities for a spectrum 
were in good agreement, so the maximum bin was easily identi- 
fied. Otherwise the radial velocities were highly dispersed, and 
this would be reflected in the cryrad of the final selected value. 

The errors on the gaussian fit of the CCF and the cr vrad values 
were used later in the analysis pipeline to identify problematic 
spectra and select alternate radial velocities where necessary, and 
also as quality criteria for the final set of results to be delivered 
to ESO. 

The S-^N jAllende Prieto et al.ll200l and lCrifoetaT](l20Tol) 
libraries provided good samples of stars with which to validate 
the radial velocities determined in the AMBRE:FEROS pipeline. 
The S"*N library comprises of 1 18 F and G dwarf stars for which 
a detailed analysis was carried out using high resolution high 
SNR spectra. All four parameters {T^s, log g, [Fe/H] and [ff/Fe]) 
as well as r adial velocities are available for each of the S'*N stars. 
ICrifo et al.l (120101) is the preliminary list of 1420 candidate stan- 
dard stars to be used to calibrate the radial velocities obtained 
with Gaia RVS. Section|6]gives further discussion on the use of 
these Ubraries to validate the AMBRE:FEROS results. 

Within the FEROS archived dataset 30 stars (338 spectra) 
were found that are also within the S'^N librar y, and 183 star s 
(411 spectra) were found that are also within Crifo et al. (2010). 
Figure and b compare the S"^N and Crifo e t al. (2010) radial 
velocities values with those determined for AMBREiFEROS. 
The mean and standard deviation of the differences between the 
two sets for each sample are also shown. There is a small offset in 
both sets wit h the S^N having a greater offset of 1.8±1.4 kms"' 
compared to .Crifo et al.i (2010) of 0.52+0.35 kms ' . The FEROS 
spectra were sampled at either 0.03 A or 0.06 A from the 
ESO:FEROS reduction pipeline coiTesponding to expected V rad 
accuracies of 1.6 kms"' and 3.3 kms"' respectively. Hence these 
offsets indic ate very good ag reement for both samples, particu- 
larly for the iCrifoetal.l(l2010i) values. 
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Fig. 7. Comparison of AMBRE:FEROS r adial velocit i es: a) S'*N 
library for 29 stars (338 spectra); and b) ICrifo et al.l (1201 Ob for 
158 stars (318 spectra). 



All of the FEROS spectra were analysed for their radial ve- 
locity. However for 42 17 spectra (19.6%) poorly defined or non- 
standard cross-coiTelation functions lead to unreliable estimates 
of the radial velocity. This is most likely due to some peculiar 
or non-stellar properties of the spectra (i.e. nova). Of these 4217 
spectra, 549 also failed the conditions of the rejection flags as 
defined in Section l4n 

Figure [8^ is a histogram of the radial velocities calculated 
for the FEROS spectra defined as good quality with well-defined 
CCF (15513 spectra). The majority of the spectra have radial ve- 
locities between -100 and +50 kms"'. Figure |8j) is a histogram 
of the FWHM of the CCF, calculated from the full CCF pro- 
file, for the same sample of spectra. It shows that the majoiity of 
spectra returned a FWHM of less than 50 kms The effects of 
the CCF on the selection of the final dataset will be discussed in 
Section |8] 

4.3. Spectral Processing B: Spectra reduction, convolution 
and first parameter estimation 

Spectral Processing B (SPB) proceeds under two iterations 
through the normalisation procedure and the MATISSE analy- 
sis. These are necessary as key tests on the normalisation and 
radial velocity correction are carried out between the two itera- 
tions allowing alterations to be made that provide better first es- 
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Fig. 8. Histogram of the radial velocity values calculated for all 
of the FEROS archived spectra defined as good spectra with 
well-defined CCFs: a) V,„<^, and b) FWHM of the Vrad CCF 



timates of the stellar parameters for the remaining pipeline pro- 
cedures. For both iterations the original archived spectra are re- 
analysed using the same procedures as outlined in SPA in terms 
of wavelength selection, cosmic ray cleaning and normalisation. 
However in both iterations of SPB the measured Nrad correction 
is automatically applied to each spectrum prior to normalisation. 
This allows each spectrum to be shifted to the laboratory rest 
frame at which the synthetic grid, and so the Bg{A) vector func- 
tions, have been calculated. 

The AMBRE:FEROS synthetic grid was set to a resolution 
lower than the resolution of the observed FEROS spectra. In or- 
der to convolve the observed spectra to the appropriate resolution 
three factors needed to be taken into account. First, the observed 
spectra are observed with constant resolving power, which cor- 
responds to increasing full-width-at-half-maximum (FWHM) of 
the spectral features with wavelength. Second, the synthetic 
spectra grid was convolved using a gaussian profile with con- 
stant FWHM hence the synthetic spectra have constant spectral 
FWHM with wavelength, not a constant resolution. Third, the 
construction of the full AMBRE synthetic spectra grid did not 
in clude astrophysical bro adening such as those due to Vsini and 
^ (Ide Lavernv et al.f2012 ^. Hence a straight-forward convolution 
using the nominal instrument FWHM and the FWHM applied to 
obtain the AMBRE:FEROS synthetic spectra grid would result 
in convolved observed spectra that had increasing FWHM with 
wavelength in disagreement with the synthetic spectra. Also for 
stars with astrophysical broadening greater than the broadening 



due to the resolving power of the instrument the spectra would 
be too broad for the appropriate parameter range of synthetic 
spectra. 

Therefore a procedure was implemented that measured the 
spectral FWHM for each of the 17 selected wavelength regions 
(see Table |2]) for each spectrum and then applied an appropri- 
ate gaussian profile to smooth each wavelength region in order 
that the convolution resulted in observed spectra that had spec- 
tral FWHM as close as possible to the constant FWHM of the 
synthetic spectra grid. 

The FWHM values were measured as part of SPB and this 
was found to significantly increase the processing time of this 
stage of the pipeline. Hence, due to it being primarily a testing 
phase, for the first iteration of SPB the observed spectra are con- 
volved using predetermined 'mean' spectral FWHM values that 
were determined for each of the 17 wavelength regions. These 
default spectral FWHM values were determined by a statistical 
analysis of a subsample of 384 of the FEROS archived spectra, 
whereby the FWHM for as many weak and medium strength 
spectral lines as possible were measured in each of the sample 
spectra. 

The FWHM was measured by fitting a gaussian function to 
each spectral feature and calculating the width of the gaussian 
profile at the midpoint of the depth of the profile. For each spec- 
trum the FWHM values increased with wavelength hence as- 
suming one FWHM for the entire spectrum was not reasonable 
across the wavelength domain of the AMBRE:FEROS analysis. 
Hence a mean FWHM value for each of the 17 wavelength sec- 
tions was determined and these were used in the convolution of 
the corresponding wavelength section. The convolution was car- 
ried out using a transformation gaussian profile for which the 
standard deviation (crx) was defined as: 



FWHMgrid^ FWHMohs^ 

0-x=A 

\ Cr fWHM Cr fWHM 

(TFWHM = V2.01n(2.0) (4) 

where FWHMofc^ is the default FWHM for the region, 
FWHMgnrf - 0.33 mA and o-fWHM transforms the FWHMs to 
crs. Therefore each region of the observed spectrum was con- 
volved such that the FWHM of the convolved spectrum matched 
the FWHM of the synthetic spectrum. 

Part of the convolution process is to resample the observed 
spectra to the same wavelength bins as the FEROS synthetic 
spectra grid. Also for SPB, the same quality tests outlined in 
SPA are used but again no spectra were rejected based on these 
tests in this first iteration. 

The convolved FEROS spectra were then analysed using 
MATISSE to obtain the first estimate of the stellar atmospheric 
parameters (T^eff, logg, [M/H] and [a/Fe]). MATISSE also out- 
puts synthetic spectra interpolated to the derived stellar param- 
eters. Within the MATISSE algorithm the determination of the 
stellar parameters undergoes 10 iterations during which the al- 
gorithm converges on the final parameters for each spectrum. 
This number is set empirically and the majority of spectra con- 
verge to their parameters in much less than 10 iterations. After 
the MATISSE analysis is complete each spectrum is tested for 
whether convergence on the final Bg(A) function occurred within 
the 10 iterations and also the spectrum is tested as to the good- 
ness of the fit between the normalised and the interpolated 
synthetic spectra in a \og(x^) calculation. In the case of non- 
convergence within MATISSE, the parameters determined at the 
ninth and tenth iterations are compared and the solution with the 
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lowest log(x^) is selected as the final set of parameters. If there 
is no convergence and/or a high \og(x^), the spectra are flagged 
for further investigation. The majority of these instances were 
attributable to incorrect selection of the radial velocity. A sepa- 
rate routine was developed that carries out an automated visual 
inspection of these spectra in order to select, where possible, an 
improved radial velocity correction from the full list of radial ve- 
locity masks. In some instances the normalisation of the spectra 
was poor due to noisy spectra, strange spectral features or an in- 
valid radial velocity correction. In the latter case the adjustment 
to the radial velocity would correct the poor normalisation. In 
the former cases these spectra were captured by quality /rejection 
flags in the normalisation process or in the construction of the fi- 
nal ESO dataset. 

After these adjustments were made the spectra were repro- 
cessed in the second iteration of SPB. This process is exactly the 
same as the first iteration except the mean spectral FWHM of the 
absorption lines for each wavelength section is now measured di- 
rectly for each individual spectrum. This second iteration of SPB 
provides the first scientific estimate of the parameters, and hence 
convolving the spectrum with spectral FWHM values measured 
specifically for each spectrum provides a more robust analy- 
sis. The measured spectral FWHM values were used to con- 
volve each region of each spectrum to match the FWHM of the 
AMBRE:FEROS synthetic spectra grid. Therefore in Equation[3] 
the FWHMofe, is the measured spectral FWHM for each wave- 
length section for each spectrum. 

Extensive testing was carried out which compared spectral 
lines that were in common between the observed and synthetic 
spectra both before and after convolution in order to ensure 
that the transformation calculation resulted in appropriately con- 
volved observed spectra. This included testing the procedure 
on the original high resolution synthetic spectra to ensure the 
AMBRE:FEROS resolution was obtained after convolution by 
this procedure. The FWHMg„v/ (0.33 mA) was calculated based 
on a sample of the AMBRE:FEROS synthetic spectra for which 
the FWHM values were determined by the same method as for 
the observed spectrum and were found to be constant with wave- 
length as expected due to how the synthetic spectra were gener- 
ated. This ensured consistency in the calculation and application 
of the FWHM in the convolution process. 

Part of the calculation of the spectral FWHM was to clas- 
sify each spectrum using the number of spectral lines and the 
widths that were measured. Three categories were used: 'Weak' 
for less than 25 lines measured which were typically found to 
have FWHM < O.Il {FWHM»,eak) and were generally noisy 
spectra; 'Medium' for more than 25 lines measured, where the 
line depths of the measured lines were between 0.5 and 0.95 
in normalised intensity, and no more than 1 strong line was 
identified (FWHM,,,^^,,,,,,); and 'Strong' for more than I strong 
line identified where the line depth of the 'Strong' line was be- 
tween 0.35 and 0.5 of the normalised intensity (FWHM„r„„j). 
The spectral lines were identified by locating the minima and 
maxima in the spectra, fitting a gaussian at each minima and 
then discarding those lines which failed the quality tests of the 
gaussian fit. For a spectrum typically classified as FWHM„,e£/,„„, 
1000 to 2000 lines were identified which were reduced to ~400 
lines with good gaussian fits. 

The 'Strong' classification was an attempt to identify the 
spectra with large, broad spectral features as this provided an- 
other estimate of the spectral classification of the star This was 
accomplished by applying a high FWHM smoothing function to 
each spectrum and then testing for which spectral features still 
remained. A spectrum without 'Strong' lines would be reduced 



to a line near the continuum, while 'Strong' features would oth- 
erwise still be prominent. This method had to be tailored to the 
FEROS wavelength region in order to place relevant limits on the 
possible number of 'Strong' lines and their location within the 
wavelength regions. If, by this method, more than 1 'Strong' line 
were found then the spectrum was classifed as FWRMsmmg and 
the mean spectral FWHM of the 'Strong' lines was recorded. 

In the cases where the spectrum was designated as 
FWHM,, 

eak or FWHMj^ong, then default spectral FWHM values 
were used instead. The spectral FWHM values for each spec- 
trum that were used in the convolution, measured or mean val- 
ues, were saved to an external file to be used in the next stage 
of the pipeline. Extensive testing was canied out using corre- 
sponding spectral lines in the original, convolved and synthetic 
spectra to confirm that the spectral FWHM procedure correctly 
convolved the archived spectra to the resolution of the synthetic 
grid across the wavelength range. 

At this stage the spectra which failed the tests identifying un- 
analysable spectra (extreme emission features, extreme noise- 
dominated spectra etc) were rejected from the analysis process. 
This reduced list was analysed once more in MATISSE in order 
to determine the first estimate of the stellar parameters to be used 
in the next stage of the analysis pipeline. The convergence and 
\ogx^ were tested again to catch any further mis-identified radial 
velocity corrections. At the end of the second iteration of SPB 
2370 FEROS spectra were rejected, 1 1% of the total number of 
FEROS archived spectra. 

4.4. Spectral Processing C: Iterative spectra normalisation & 
parameterisation 

The final number of FEROS spectra that were passed to the final 
stage of the analysis pipeline. Spectral Processing C (SPC), was 
19181. SPC again uses the same procedures as SPA and SPB 
with two differences. The first difference is that the spectra are 
convolved using the spectral FWHM values determined in the 
second iteration of SPB rather than re-measuring the FWHM. 
The second difference is the use of a more robust method for the 
normalisation of the spectra. 

Previously, in SPA and SPB, a 'rough' normalisation pro- 
cess was carried out that used pre-selected continuum regions to 
normalise the spectra to unity. But this is an invalid assumption 
for many stars, in particular for cool stars which have depressed 
continuum regions due to molecular bandheads that should not 
be normalised to unity. However at the beginning of SPC we 
have an estimate of the stellar parameters from the second iter- 
ation of SPB that gives us some valid information about each 
spectra that we did not previously possess. 

We take advantage of this information by using the interpo- 
lated synthetic spectra generated for each set of stellar param- 
eters as a better estimate of the continuum placement for the 
normalisation of the observed spectra. In an iterative process 
that discards absorption and emission features by polynomial fit- 
ting and sigma clipping to leave only 'continuum' regions, the 
synthetic spectrum is divided out of the corresponding observed 
spectrum leaving a residual of continuum points. A continuum 
profile is fitted to this residual which is then divided out of the 
observed spectrum, normalising it to a pseudo-continuum that 
better represents the parameters of the star. This process over- 
rides any residual curvature that remains after, or was introduced 
by, the normalisation to the hot star continuum. 

These pseudo-continuum normalised spectra are reanalysed 
in MATISSE to obtain a second estimate of the stellar parameters 
and another set of synthetic spectra. The 'normalisation treat- 
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Fig. 9. Normalisation procedure for a two stars with T^ff ~ 4600 K (LEFT) and T^fi ~ 5300 K (RIGHT) shown by the wave- 
length region about two lines of the magnesium triplet. TOP: the second iteration of SPB, MIDDLE: first iteration of SPC, 
and BOTTOM: final solution from SPC. The stellar parameters (T^ff, \ogg, [M/H] and [ff/Fe]) and goodness of fit over all the 
AMBRE:FEROS wavelengths (log;^'^) are shown. The synthetic spectrum interpolated to these parameters is shown in red, while 
the normalised observed spectrum is shown in black. The difference between observed and synthetic spectra is also shown at 0.0. 



ment & MATISSE analysis' cycle is repeated 9 times, which had 
been determined to be a sufficient number of iterations within 
which the normalised spectra and stellar parameters could con- 
verge to their optimum state. This typically occurred within five 
iterations. This approach provides a robust incremental adjust- 
ment of the parameters and normalisation process that hones in 
on a realistic estimate of the stellar parameters of each spectra, a 
process which is not possible using the normalisation method of 
SPA and SPB. 

Figure |9] shows the normalisation procedure for two FEROS 
spectra in the wavelength region about two of the magnesium 
triplet at 5160A as an example of depressed continuum regions. 
The two stars correspond to temperatures of Teff ~ 4600 K 
(SNR = 115) and T^ff ~ 5300 K (SNR = 204) respectively. For 
each spectrum the normalisation by SPB, the first iteration of 
SPC and the final solution that was found in SPC are shown. 
The stellar parameters and \ogx^ (which is calculated over all 
the AMBRE:FEROS wavelengths) determined at each stage are 
shown as well as the synthetic spectrum interpolated to the stel- 



lar parameters. The difference between the observed and syn- 
thetic spectra is also included to show how it converges. 

For the cooler star this rough normalisation of SPB is not 
a good fit as the observed spectrum is set too high compared 
to the synthetic spectrum generated at the corresponding stellar 
parameters. This shows how inadequate normalisation to unity 
is for cool star spectra. However the rough normalisation does 
provides a reasonable fit to the synthetic spectrum for the hot 
star although the placement of the continuum is still not ideal. 

At the initial iteration of SPC the spectrum undergoes its first 
normalisation to the previous solution synthetic spectrum. For 
the cooler star the observed spectrum now sits below unity in 
better agreement with the new synthetic spectrum, but there are 
still mismatches in the relative placement of the two spectra and 
the line depths. For the hotter star there is very good agreement 
between the observed and synthetic spectra. In both cases the 
improvement in the match are shown by the lower log;^'^ values. 

The final solution for the cooler star (/ = 6) shows excel- 
lent agreement between the observed and the synthetic spectra, 
which is reflected in the much lower \ogx^ value. Indeed the vi- 
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sual inspection of the spectra show that they are in good agree- 
ment in this section of wavelength. For the hotter star the solu- 
tion (/ = 7) is very close to the solution determined in the first 
iteration of SPC, indeed the log;^'- does not change. Visually this 
good fit is obvious between the observed and synthetic spec- 
tra. The comparison of these two stars illustrates the greater 
difficulty there is in the normalisation and parameterisation of 
cool stars. However the procedure developed for the AMBRE 
pipeline can successfully manage both cases and the normalisa- 
tion procedure in SPC is particularly effective in the case of cool 
stars where the continuum regions are more heavily obscured by 
spectral features. 

Typically convergence occured by the fourth or fifth iteration 
of the normalisation and stellar parameter determination cycle. 
When convergence has occurred the final set of parameters are 
selected as those from the converged solutions with the lowest 
logCf^). In the case of non-convergence within the 9 iterations 
of SPC, it is assumed that for the final six iterations any bias 
introduced by the initial rough normalisation has been erased. 
Of these final six iterations the solution with the lowest log{x^) 
is selected as the final stellar parameters for the spectra along 
with the corresponding normalised and synthetic spectra. 

Hence the radial velocity, stellar parameters (T^ff, \ogg, 
[M/H] and [c/Fe]), normalised spectra and corresponding in- 
teprolated synthetic spectra were determined in SPC for 19181 
of the 21551 FEROS spectra. 

Extensive testing of the FEROS analysis pipeline was car- 
ried out in order to optimise the procedures and to validate the 
MATISSE results by making comparison to literature values. 
The following sections describe the different tests and valida- 
tions that were carried out. 



5. AMBREiFEROS internal error analysis 

The internal errors can be used to determine how well MATISSE 
will derive the stellar parameters of similar stellar types which 
have different levels of noise and different uncertainties in radial 
velocity. In order to test for noise and radial velocity effects a 
sample of 500 synthetic spectra were generated at random stel- 
lar parameters covering the entire stellar parameter range, spec- 
tral resolution and wavelength domains of the AMBRE:FEROS 
synthetic spectra grid. 



5.1. Internal errors: SNR 

To test the effects of noise the sample of synthetic spectra was re- 
produced five times by adding differing levels of noise per pixel 
(0.1 A per pixel for AMBRE:FEROS synthetic grid) at SNR = 
10, 20, 50, 100. Each sample was then re-analysed in MATISSE 
and the derived stellar parameters were compared to the original 
values. 

Figures [TOh . to d. show how the difference in stellar param- 
eters changes with increasing SNR for Teif, \ogg, [M/H] and 
[or/Fe] respectively. For each noise-added sample the difference 
between the original and derived parameters was calculated and 
then the difference value at the 70th percentile (~ Icr) was deter- 
mined. For example, in FigurefTOk. at a SNR of 10 (the noisiest), 
70% of the synthetic spectra returned Teff values within 14 K of 
the original T^g value. At a SNR of 100, 70% of the sample 
returned values within 2 K of the original T^ff value. Similarly 
small errors were also found for log g, [M/H] and [a/Fe] that di- 
minished to negligible levels at SNR of 100 in all cases. This 
analysis shows that for two synthetic spectra with very similar 
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Fig. 10. Internal error for each parameter with changes in SNR. 
a) The Areff that 70% of the synthetic sample were less than or 
equal to when calculating the difference between the nominal 
Teff and the Teff determined for the respective SNR. b) As for a) 
but for log g. c) As for a) but for [M/H]. d) As for a) but for 
[a/Fe]. 



intrinsic stellar parameters, where the differences in the spectra 
are only attributable to noise, MATISSE will determine close 
to the same stellar parameters indicating that the internal errors 
of MATISSE are very small. This can be attributed to the large 
wavelength coverage and high resolution of the FEROS spec- 
tra which contain a great deal of information on the stellar pa- 
rameters in terms of spectral lines, even at low SNR, of which 
MATISSE takes advantage. 

Overall this analysis shows that the internal error due to SNR 
is negligible in the FEROS results. The internal (int) error on 
each parameter due to SNR {(Tint^snr) is calculated as part of the 
MATISSSE program. It was used to calculate the total internal 
error as explained at the end of this section. 

5.2. Internal errors: Vrad &. Vsini 

A similar investigation was carried out in order to determine 
the effect of uncertainties in radial velocity on the determina- 
tion of the stellar parameters by MATISSE. The sample of in- 
terpolated spectra were reproduced this time undergoing radial 
velocity shifts of AV^rf - 0, 0.5, 1.0, 5 and 10 kms"' with an 
SNR of 100. Figures [TT^ to d show how the difference in stel- 
lar parameters changes with increasing radial velocity uncer- 
tainty for Teff, \ogg, [M/H] and [a/Fe] respectively. Again the 
difference between the original and derived stellar parameters 
was obtained and the difference value at the 70''' percentile for 
each sample was calculated. At a radial velocity uncertainty of 
5 kms"' (~ 0.09 pixels) the error in stellar parameter determi- 
nation becomes significant (>1%) for each of the stellar param- 
eters. However if the uncertainty in the radial velocity is less 
than 5 kms ' the MATISSE parameters compare well to the true 
values. 
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the instrument. The AMBRE synthetic spectra grid was gener- 
ated with no variations in Vsini, assuming all stars to be slow 
rotators. However it is important to consider whether variations 
in Vsini will have an effe ct on the determined stellar parame- 
ters . iGazzano et al. carried out a robust investigation into 
the Vsini limits of reliable stellar parameter determination using 
MATISSE and found that spectra with Vsini < 11 kms ' pro- 
duced good results in th e MATISSE analysis o f a sample of FGK 
dwarfs. Based on this, iGazzano et al.l (1201 Ol) accepted all CCF 
FWHM < 20 kms^^ for spectra at a resolution of R~26,000. 
The FEROS sample and the AMBRE synthetic spectra grid are 
of much broader range in stellar parameter s and waveleng t h, and 
of lower resolution (R~15,000) than the Gazzano et alJ (l2010h 
study. Hence the effects of Vsini on the range of stars for which 
the synthetic grid has been optimised, would be well-masked for 
the AMBRE:FEROS configuration. In Section [7.2.3l we examine 
the FEROS sample in the context of the CCF FWHM to identify 
spectra which are not well-represented by the AMBRE synthetic 
spectra grid. 



5.3. Internal errors: Normalisation 



Fig. 11. Internal error for each parameter with changes in the 
Vrad uncertainty (kms"'). a) The Areff that 70% of the synthetic 
sample were less than or equal to when calculating the difference 
between the nominal Teff and the Teff determined for the respec- 
tive Wrad uncertainty, b) As for a) but for log g. c) As for a) but 
for [M/H]. d) As for a) but for [a/Fe]. 

The internal uncertainty in each stellar parameter due to ra- 
dial velocity was calculated using the equations that connect the 
points in Figure [TT] The error on the radial velocity, crCVrad), 
determined in the radial velocity programme was used as the in- 
put to determine the internal error due to the V rad uncertainty 
(o'int.Vrad) foi" each parameter for each spectrum according to the 
following conditions: 



The effects of normalisation on the MATISSE determination 
were measured using the iterative normalisation process (SPC) 
of the FEROS analysis pipeline. Using a test sample of 384 
FEROS spectra ((SNR)=100+45) the changes in normalisation 
between iterations of SPC was explored. 

Figure [12] shows the progression of the goodness of fit be- 
tween the reconstructed spectra and the normalised spectra at 
each iteration of SPC as a Box & Whisker graph. Assuming that 
by the 9''' iteration convergence has occured, the absolute differ- 
ence in the \ogx^ between the 9''' and i''' iteration was calculated 
for each of the 384 test spectra. The spread in values is largest at 
/ = and there is a large decrease in the spread for / = 1 . The 
following iterations show a more sedate decrease in spread and 
appear to have converged from ; = 6. 
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The rotational velocity {Vsini) of a star can have an impact 
the broadening of the observed spectral features depending on 
the magnitude of the rotational velocity and the resolution of 
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Fig. 12. Box & Whisker graph of the absolute difference in the 
log;^^^ of the reconstructed to the normalised spectra for the 384 
test spectra between the 9''' SPC iteration and each of the previ- 
ous iterations (| \ogx^(9) - log;^'-(OI). The box is constrained by 
the 25''' and 75''' percentiles, and the red lines are the medians. 
The whiskers extend to the furthest data points not considered to 
be outliers. 
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Fig. 13. Internal error with changes in Normalisation for a) Teft; b) \ogg; c) [M/H]; and d) [a/Fe]. The red lines are determined 



as the 2cr uncertainty on values in two SNR bins (ASNR 
125. 



30) centred at SNR = 65 & 100 and constant uncertainty for SNR > 



To measure the internal error associated with the normalisa- 
tion process we investigated the difference in stellar parameters 
between the 9''' and 7''' iteration for the test sample. Hence we 
assumed that the solutions have converged from at least the 7''' 
iteration and the subsequent variations in the stellar parameters 
are due to minor adjustment of the normalisation of the observed 
spectra. The variation from convergence (/ = 7) to the final it- 
eration (/ = 9) was used to give the internal error due to the 
normalisation process. 

Figure [13] shows the difference in each stellar parameter be- 
tween the 9''' and 7''' iterations against SNR for the 384 test spec- 
tra as black points. These values were binned into two SNR bins 
(ASNR « 30) centred at SNR = 65 & 100, and a single bin for 
all points with SNR > 125. The 2 cr uncertainty for each bin is 
shown as red squares at the bin centre. Red lines connect the bin 
uncertainty values. For T^ff the spread is constant for each bin 
and so a constant value for the internal uncertainty due to nor- 
malisation was adopted. For log g, [M/H] and [a/Fe] the spread 
diminishes with increased SNR until SNR =125 after which a 
constant spread is found. The adopted internal errors for these 
three parameters were therefore a piecewise function of SNR. 
The following equations define the internal errors due to nor- 
maUsation for each parameter (cr(0),„, „o„,i): 



ll.OK 



-0.0017 SNR + 0.2414 dex if SNR < 100, 
o-(log ^;),„,.„o™ = I -0.0018 SNR + 0.2500 dex if 100 < SNR < 125, 
-0.025 dex if SNR > 125, 



o-m/H]),,, 



o-ila/Fe]),,, 



-0.0007 SNR + 0.0964 dex 
-0.0004 SNR + 0.0690 dex 
-0.014 dex 

-0.0006 SNR + 0.0671 dex 
-0.0002 SNR + 0.0260 dex 
-0.006 dex 



if SNR < 100, 

if 100 < SNR < 125, 

if SNR > 125, 

if SNR < 100, 

if 100 < SNR < 125, 

if SNR > 125, 



5.4. Total internal error 

The internal errors due to the SNR, V rad and NormaUsation cal- 
culated above were combined in quadrature to give the total in- 
ternal error for each parameter (cr(0)„„) for each spectrum. This 
was calculated as follows: 



Cr(0)i„, = ^Cr'^{9)inl,snr + CT^O) i„,,y + Cr^{,9)int,norm (5) 

These values were reported for each spectrum in the ESO 
dataset as described in Section[8] 



6. AMBREiFEROS external error analysis 

The external error was quantified by comparing the 
AMBRE:FEROS stellar parameter values to literature val- 
ues for key reference stars that exist within the AMBRE:FEROS 
dataset. 
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A list of stars that had been observed with FEROS but were 
also part of other high quality spectroscopic studies were iden- 
tified and this reference sample was used to provide quality cri- 
teria for each derived stellar parameter. Unfortunately, reference 
stars covering the whole range of possible parameters could not 
be found in the literature. In several cases, estimates of the mean 
stellar metallicity were lacking and the situation was even worse 
for the [a/Fe] chemical contents since too few reference stars 
with known [a/Fe] content could be found in the literature. 

6.1. Building the FEROS reference sample 

In order to compare the MATISSE stellar atmospheric parame- 
ters to the literature and to define some quality criteria, a list of 
reference stars was built using key spectral atlases and libraries. 

6.1.1. Stellar Atlases: Sun, Arcturus, Procyon 

Table |4] lists the standard stars and the corresponding stellar at- 
lases that were used to test the FEROS analysis pipeline. Of these 
atlases only the S'*N solar atlas was observed using FEROS. 
However the adaptability of the pipeline meant that the use of 
spectra from other instruments presented no difficulty as long 
as the required wavelength ranges were available. The accepted 
stellar parameters for each standard star are listed in bold in 
Table |4] The final stellar parameters that have been derived in 
the AMBRE:FEROS analysis with bias cotTections (see end of 
this section) for each atlas are also given. 

Table 4. List of standard stars and the corresponding spectral 
atlases used to test the FEROS analysis pipeline. The accepted 
parameters for the standard stars are listed (bold) as well as the 
final stellar parameters with bias corrections derived for each 
atlas in the AMBRE:FEROS analysis. 



Standard star/ATLAS 




^ogg 


[Fe/H] 


[a/Fe] 




(K) 


(dex) 


(dex) 


(dex) 


SUN 


5770 


4.44 


0.00 


0.00 


Wallace et al. (1998)" 


5732 


4.42 


-0.03 


0.05 


Brault & Neckel (1999)'' 


5674 


4.32 


0.00 


0.08 


BASS2000'' 


5841 


4.49 


-0.01 


0.06 


S-'N'' 


5735 


4.43 


-0.03 


0.05 


ARCTURUS* 


4300 


1.70 


-0.60 


0.20 


Hinkle et al. (2000r 


4537 


2.12 


-0.49 


0.28 


UVES-POP" 


4468 


2.03 


-0.47 


0.29 


PROCYON* 


6600 


4.00 


0.00 


0.00 


UVES-POP'' 


6538 


4.00 


-0.25 


0.07 



ftp://ftp.noao.edu/catalogs/arctumsatlas/visual/ 
ftp://ftp.hs.uni-hamburg.de/pub/outgoing/FTS-Atlas/ 

''■'Delbouille et al. (1973); http://bass2000.obspm.fr/ 

AUende Prieto et al. (2004); http://hebe.as.utexas.ed u/s4n/| 
Bagnulo et al. (2003) ; http ://w w w. sc .eso . org/ san tiago/uvespop/ 1 

*** Mean of parameters listed on SIMBAD 

These atlases were combined in a test sample that in- 
cluded high resolution synthetic spectra for the Sun and 
ArctutTis, as well as synthetic spectra for each grid point in the 
AMBRE:FEROS synthetic grid corresponding to the Sun and 
Arcturus. These were used to check the internal consistency of 
the pipeline in terms of normalisation and convolution of the 
input spectra. The results for the Sun are in very good agree- 



ment with the accepted values. For Arcturus we over-estimate 
the Tgff and log g values however we are within 2 cr of the ac- 
cepted values for Teff. For Procyon there is good agreement 
in Teff and \ogg, although the T^ff is slightly under-estimated 
within the stated external errors. However the metallicity is 
under-estimated by -0.25 dex which is outside the 2 cr limit. 
The under-estimation of the T^ff may account for some of this 
discrepancy in metallicity. However further investigation into 
understanding the variation in the stellar parameters for these 
standard stars is ongoing. 

6.1.2. Spectral libraries reference sample 

Preliminary investigations showed that merely making compar- 
ison between the MATISSE results and all of the stellar param- 
eter values that were available in SIMBAD for the FEROS stars 
resulted in a great deal of disparity. There was no quality de- 
termination on the values with which we made the comparison. 
Homogeneously analysed large sample datasets were necessary 
in order to ensure that a valid comparison was being made. To 
this end key spectral libraries were investigated in order to locate 
within the FEROS archive dataset stars that had been analysed 
in high quality studies. This list of spectral libraries is given in 
Table|5] 

The p rimary spectral library t hat was used was the PASTEL 
database (ISoubiran et al.ll201 ()'FL which provides a selection of 
key papers that carried out detailed spectroscopic studies at high 
resolution and high SNR to derive high quality stellar parame- 
ters. From these papers we were able to identify a total of 148 
stars cotTesponding to 618 FEROS archived spectra (Table |5). 
There was significant crossover of stars between the studies 
which provided another level of comparison for this work. From 
the PASTEL sample ~120 stars (~300 spectra) had values for 
three stellar parameters, T^ff, log g and [Fe/H], while the remain- 
ing stars had T^ff values only. In the following comparison we 
assume the literature values to be correct. However we note that 
the literature values, while they are high q uality stud ie s, wer e 
determined using a range of techniques (fSoubiran et aDl2010l) . 
in particular, several studies determined photometric T^ff val- 
ues using the Infra-Red Flux Method ( Alonso et al. 1996, 199^^ 
'Ran urez & Melended l2005t iGonzalez Hernandez & Bonifacia 
200^ . 

Figure [T4h. to b compare the T^ff, logg, [M/H] and [a/Fe] 
values determined by AMBRE:FEROS with the literature val- 
ues of the reference sample and stellar atlases. The legend in 
each provides the bias and cr for each subsample. The axes rep- 
resent the limits of the final accepted parameters and clearly 
the reference sample does not cover the entire grid as would be 
ideal. Despite this limitation the reference sample was sufficient 
to define bias corrections which were necessary to apply to the 
AMBRE:FEROS stellar parameters. 

6.1.3. Reference sample: Teff 

Figure [T4h . compares the Teff for each star It is important to 
note that many of the literature T^ff values were not determined 
spectroscopically, as stated above. Despite the differences in 
measurement technique, overall there is excellent agreement be- 
tween the AMBRE:FEROS results and the PASTEL values, and 
also good agreement for the stellar atlases. There is a slight turn- 
ing over of the distribution from T^ff >~ 5750 K for which a 



^ http://pastel.obs.u-bordeaux 1 .fr 
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Table 5. List of spectral libraries used to select comparison stars from the FEROS archived spectra. The number of stars from each library and the corresponding number of 
FEROS spectra are also listed. 



Source 



Gaia RVS Standards 



PASTEL 
Soubiran et al. 



Papers 



Abbrv. 



No. Stars No. Spectra Methods applied 



Allende Prieto et al. (2004) (reff,logg,[Fe/H],[a/Fe],V,.„rf) 

Crifo et al. (2010) (Y.ad) 

BensbxetaL Q003) (Jeff Jog g, [Fe/H] , [a/Fe] )" 

reff,logg,[Fe/H] 

Fuhrmann et al. (1997, 1998); Fuhrmann (1998a,b, 2004, 2008) 
Gratton et al. (1982); Gratton (1989); Gratton et al. (1996, 2003) 
Hekker & Melendez (2007) 

Luck & Bond (1983, 1991); Luck & Heiter (2005. 2006) 

Mc William (1990) 

Mishenina & Kovtvukh (2001) 

Mishenina et al. (2003, 2004 . 2006. 2008) 

Ramirez et al. (2007)" 

Valenti & F ischer (2005) 

Tea only 

Alonso et al. (1996, 1999) 
Blackwell & Lynas-Gray (1998) 
di Benedetto (1998) 

^onzalez Hernandez & Bonifacio (2009) 
Kovtvukh et al. (2003. 2004 . 2006) 
Masana et al. (2006) ' 
Ramirez & Mele ndez (2005) 

PASTEL Total 



AP04 



CIO 



30 



158 



338 



318 



Alonso et al. (1996, 1999) calibration, isochrones, 
Allende Prieto (2003) GA, spectroscopic inversion 
Convolution, CORAVEL cross-correlation 



FEROS Observations 2000-2001 
GA = Genetic Algorithm, EP = Excitation Potential, IP = lonisation Potential 



B03 


68 


68 


Fe I EP balance, parallax, equivalent widths 


F97-08 


49 


206 


Balmer line wings, Fe IP balance, Fe I/Fe II line profiles 


G82-03 


23 


118 


Fe I EP balance,Fe I/Fe II IP balance, equivalent widths 


H&M07 


6 


27 


Fe I EP balance, Fe I/Fe II IP balance, equivalent widths 


L&H83-06 


41 


224 


Fe I EP balance, Fe I/Fe II IP balance, equivalent widths 


McW90 


10 


41 


IRFM+colour calibration, parallax, equivalent widths 


M&K01;M04-08 


42 


227 


Wings of H^, Fe I EP balance. 








Fe I/Fe II IP balance, equivalent widths 


R07 


26 


151 


Colour calibration, parallax, equivalent widths 


V&F05 


81 


310 


'Spectroscopy Made Easy' 


A96,99 


22 


147 


Infra-Red Flux Method (IRFM) 


B&LG98 


20 


152 


IRFM 


dB98 


18 


163 


Surface brightness and V - K calibration 


GH&B09 


26 


162 


Line depth ratios 


K03-06 


48 


354 


SED fit + V - A" calibration 


M06 


57 


306 


IRFM 


R&M05 


37 


197 


IRFM + colour calibration 




148 


618 
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Fig. 14. Comparison of stellar parameters determined in AMBRE:FEROS (A:F) and the reference sample (RS) from PASTEL: a) 
Tgff, b) log^, c) [M/H] compared with [Fe/H], d) [a/Fe]. The legend identifies the relevant individual papers and the associated 
sample bias+cr. The stellar atlases are also shown. The dashed line indicates a 1:1 agreement. The axis limits represent the limits 
on the final accepted parameters. 



small bias correction was applied that will be quantified at the 
end of this section. 



tween [Fe/H] and [M/H] then made a systematic bias correction 
to the final [M/H] results as outHned at the end of this section. 



6.1.4. Reference sample: log g 

Figure [T4b. shows the comparision of the \ogg values for which 
there is reasonably good agreement. However the majority of 
the sample are dwarfs (high log g), and the giants (low log g) are 
under-represented. There is a slight gradient within the sample 
which resulted in a small bias correction being applied at high 
logg for the final results (see end of the section). The lack of a 
significant sample of giants made it difficult to accurately define 
a bias at low log g, hence the correction was only applied at high 
logg. However the uncertainty on the giants is definitely higher 



6.1.5. Reference sample: [M/H] 

Figure fT4b . shows the comparison of the AMBRE:FEROS 
[M/H] values with the literature [Fe/H] values. The difference in 
definition between [M/H] and [Fe/H], as discussed previously, 
means that this is not an accurate comparison, as so many more 
elements (albeit of lesser contribution) than just Fe are included 
in the MATISSE metallicity. The comparison shows an over- 
all systematic offset between the AMBRE:FEROS and literature 
values reflected in the subsample biases. The sources of this bias 
are difficult to quantify so we assumed a direct comparison be- 



6.1 .6. Reference sample: [a/fe] 

Unfortunately very few comparison stars with published values 
of their [a/Fe] e xist in the literature. The stellar parameters in 
the S'*N library (lAUende Prieto et al.ir2004 ) also include a el- 
ement abundances where possible, which provided some com- 
parisons with the AMBRE:FEROS [a/Fe] results. The S'^N T^ff 
values were determined using photometric calibrations while the 
log g were determined from Hipparcos parallaxes. The chemical 
abundances were determined using minimisation of the spec- 
tral line profile b etween the observed spec trum and a grid of 
synthetic spectra (' Allende Prieto et al ] l200l . 

A key study that was investigated specificall y in order to 
valida te the AMBRE:FEROS [a/Fe] values was Bensbv et al.l 
(l2003h . This study was a detailed analysis of 66 F and G dwarf 
stars in the galactic disk, and accurate abundances were de- 
termined for the a elements, Mg, Si, Ca, Ti. The mean of 
these abundances was taken as the value for the global [a/Fe] 
with which to compare the AMBRE:FEROS [a/Fe] values. This 
study was particula rly useful because the spectra analysed in 
i Bensby et all (l2003l) were observed with FEROS. However the 
observations took place in 2000 and 2001 and so the spectra 
were not part of the archived sample delivered to AMBRE. The 
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original spectra were obtained (Bensby, private communication) 
and analysed in the AMBRE:FEROS pipeline. Hence a direct 
comparison between AMBRE:FEROS and Bensby et al. (200l) 
could be made for all four parameters (See Figures HAk to c). 



Figure [T4t l. shows the comparison of the AM BRE:FEROS 



i g/Fe" results with the values from the S'^N librarv. ^ensbv et al] 
2003h and the stellar atlases. The biases are also listed for 
each sample. There is reasonably good agreement with these 
samples. In particular the r e is e xcellent agreement in the re- 
sults for the iBensbv et al.l (l2003l) stars, which as stated above, 
is a study where the a element abundances were carefully 
determined, thereby providing an excellent validation of the 
AMBRE:FEROS [a/Fe] results. 



6.3. External error 

The external error for each stellar parameter was determined 
from the above reference sample analysis. For each individual 
study the mean difference and spread in differences (cr) between 
the literature and AMBRE:FEROS values was calculated. The 
mean of these cr values for each stellar parameter was taken as 
the global external error ((T(ff)ext) for all of the spectra as follows: 



cr(reff)„, - 120 K 

f0.20dex iflog?>3.2, 
[0.37 dex iflog^<3.2, 

cr([M/H])„, = 0.10dex 

cr([Q'/Fe])„, = 0.10 dex 



6.2. B/as corrections 

The sample of reference stars and the spectral atlases pro- 
vided a crucial comparison at all stages of the development 
of the AMBRE:FEROS analysis pipeline. Although it was not 
possible to cover the entire parameter space, and the refer- 
ence sample was ultimately biased towards metal-rich dwarfs 
over a small temperature range, the results provided sufficient 
information with which to identify biases within the analy- 
sis, assuming the literature values to be correct. The follow- 
ing corrections were made in order to remove these biases: 

I T^ff - 35 if Teff < 5300 K, 

T^sicor) = \T,ff + 0.21 x T^g - 1141.4 if 5300 < < 6000 K, 
[reff+UO if Teff > 6000 K. 

log g - 0.296 X logg + 1.388 if log g > 4.0, 

, , logg 0.204 if 3.75 < logg < 4.0, 

log s(cor) = ■{ 

] log g 0.8 17 X logg - 2.860 if 3.5 < logg < 3.75, 

logg iflogg<3.5. 

[M/H](cor) = [M/H] +0.15 
[a/Fe](cor) = [a/Fe] 

There are several potential sources of the bias corrections for 
the AMBRE:FEROS analysis. First, the diff'ering techniques, as 
well as differing spectral domains and resolutions, used in the 
determination of the stellar parameters by the reference stud- 
ies compared with AMBRE:FEROS analysis may make a sig- 
nificant contribution to the difference in derived parameters. 
Second, while the normalisation process within the AMBRE 
analysis pipeline was designed to be as robust as possible, the 
normalisation of stellar spectra is an inherently complex prob- 
lem, particularly over large wavelength ranges with many spec- 
tral features. For spectra of particular spectral types (e.g. cool, 
low gravity, metal-rich) the normalisation procedure may not be 
as robust as for less detailed spectra. Third, due to the length 
of the atomic and molecular linelist that was required to synthe- 
sise the AMBRE:FEROS wavelength regions, it was inefficient 
to carry out a calibration of the linelist to standard stars (i.e. the 
Sun). It was assumed that the number of lines was statistically 
sufficient to dampen the noise from ill-fitted spectral features 
but non-calibration of the line list may play some role in the de- 
gree of the bias corrections. Individually, and in combination, 
these are the most likely sources of the bias corrections for the 
AMBRE:FEROS stellar parameters. 



7. ESO Table: Rejection criteria 

The construction of the ESO Table of stellar parameters that 
would be delivered to the ESO Archive was the final stage of 
the pipeline. At this stage the quality flags and tests included 
throughout the pipeline were drawn together to provide the fi- 
nal set of stellar parameters. The columns headings, their defini- 
tions, range of values, null values and rejection conditions used 
to construct the ESO Table are listed in Table lA.fl 

There were two phases of rejecting spectra in the 
AMBRE:FEROS analysis. First, prior to analysis in SPC spec- 
tra were rejected due to spectral quality issues. This will be de- 
scribed in the Section 17.11 Second, the remaining spectra were 
all analysed in SPC for their stellar parameters and then another 
set of rejection criteria were applied to construct the final table. 
These criteria are described and discussed in Sections [7 .2 I to 17 .41 

7. 1 . Pre-SPC: Non-standard spectra 

Prior to SPC, 11% of the FEROS archived spectra were re- 
jected as being non-analysable primarily due to being non- 
standard spectra. The rejection flags, as defined in Section l4~n 
were for 'Faulty Spectra', 'Extreme Emission Features', 'Poor 
Normalisation' and 'Excessive Noise'. Often the nature of a 
spectrum rejected at this stage meant the conditions for two or 
more of these flags were met. For instance, an extremely noisy 
spectrum was also consequently poorly normalised. Quality 
flags were also attached for 'Large Emission Features' and 
'Instrumental Relics' but these spectra were not rejected based 
solely on these flags. The number of spectra identified for each 
rejection flag and for each quality flag are listed in Table |6] 

7.2. Rejection criteria post-SPC 

The remaining 89% of the FEROS spectra were analysed in SPC 
and so stellar parameters were obtained for each of these spectra. 
Figure [TSh shows the HR diagram for all the stellar parameters 
obtained in SPC. The giant branch and main sequence can be 
observed but there is a great deal of mis-classification at hot tem- 
peratures for log g > 2 dex, and an overdensity at Teft ~ 3000 K 
and log g ~ 2 dex. Most noticeable is that there are many spectra 
(8947) for which the derived parameters lie outside the rgft and 
log g grid limits which are indicated in red. In the construction 
of the ESO Table it would have been very straightforward to just 
reject all spectra whose parameters lay outside the limits of the 
synthetic grid and so only report to ESO those spectra with pa- 
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Fig. 15. HR diagrams of the FEROS stellar parameters as each rejection criterion is applied as follows: a) All spectra analysed 
in SPC; b) Spectra with M ,ad CCF with positive contrast, ^ < 0.2 and fgg- < 0.1; c) o-y,,,, < 10 kms^i; d) CCF FWHM < 

40 kms e) FWHM of medium strength spectral lines < 0.8 mA; f) FWHM of strong spectral lines < 8 mA. The Jeff and logg 
boundaries of the synthetic grid are indicated in red. 



Table 6. The number of spectra rejected (Rej) prior to SPC based 
on the defined rejection flags and the number of spectra accepted 
(Acc) for SPC but identified by the quality flags. 



No. % of % of 

Spectra FEROS Rej/Acc Sample 



FEROS 


21551 






Rejected before SPC 


2370 


11 




Faulty Files 


71 


0.33 


3.0 


Extreme Emission Features 


1455 


6.75 


61.4 


Excessive Noise 


109 


0.51 


4.6 


Poor Normalisation 


735 


3.41 


31.0 


Accepted for SPC 


19181 


89 




Large Emission Features 


3177 


14.7 


16.6 


Instrumental Relics 


58 


0.27 


0.30 



rameters inside the grid. However to have the synthetic grid lim- 
its as the sole reason to reject so many spectra (~ 47%) without 
further explanation seemed inadequate. Thus we further investi- 
gated the observational and astrophysical characteristics (indica- 
tors) of the spectra in order to better understand the dataset and 
to be able to impose more robust quality control criteria on the 
final dataset to be reported to ESO. 



The SNR and logx^ were investigated as potential indica- 
tors but both proved to be functions of spectral type (the more 
spectral features the lower the SNR, the poorer the observed to 
synthetic spectra fit) and the criteria rejected RGB stars rather 
than the high temperature mis-classifications. 



Indicators relating to the measurement of the spectral 
FWHM and the radial velocity proved to be the least affected 
by spectral type and better identified the mis-classifications. The 
key indicators that were explored were: the Wrad CCF contrast, 
the error on the amplitude and the error on continuum of the V,-,,^ 
CCF; the V™^ error; the FWHM of the Yrad CCF; the spectral 
FWHM; and the AMBRE:FEROS synthetic grid boundaries in 
Teif, logg, [M/H], [a/Fe]. 



Table |7] summarises the number of spectra that satisfied the 
conditions of each rejection criteria for the spectra rejected be- 
fore and after SPC. The post-SPC rejection criteria were applied 
to the spectra rejected before SPC to further characterise that 
sample. To illustrate the effects of the rejection criteria on the 
sample of spectra analysed in SPC, Figure [15] shows progres- 
sively the HR diagram of the AMBRE:FEROS stellar parame- 
ters as each criteria is applied. The limits of the synthetic grid 
in Teff and log g are also shown. This process is discussed in the 
following sections. 
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Table 7. Summary of the number of spectra flagged as satisfying 
the conditions of the post-SPC rejection criteria. The spectra that 
were rejected prior to the SPC analysis were rejected due to other 
quality issues (see Section 17. 11 1. The spectra analysed in SPC 
were accepted or rejected based on the listed criteria. 





Spectra rejected 
pre-SPC 


post-SPC 


Total 


2370 


12673 


PD CCF" 


549 


3668 




309 


6348 


CCF' 


469 


10062 


FWHM'' 


7 


1199 


Grid^ 




8494 



a. Poorly-Defined CCF (Se ctionlTm 

b. o-v,.„^ > 10 km/s (Section|T2j2} 

c. FWHM CCF > 40 km/s (Section lT231 l 

d. FWHM,„frf„„„ > 0.8 mA or 
FWHM.„,,,„j > 8 mA (Section|T231l 

e. Parameters lie outside grid parameter space (Section IT? 



7.2.1. CCF contrast, amplitude & continuum 



As described in Section 14^21 the radial velocity was determined 
using binary masks to calculate a CCF. The contrast, the depth 
(amplitude) of the CCF, the continuum placement of the CCF 
and their associated errors (as calculated in the IDL:GAUSSFIT 
routine) are measures of the quality of the CCF. The rejection 
criteria were derived from these quantities as follows: 

1 . A negative contrast means that the CCF profile is inverted, 
which is contrary to the expected result, hence all spectra 
with a negative contrast were rejected. 

2. The relative error on the amplitude of the CCF ij^) is 
mainly a measure of how much noise there is in the CCF. If 
the < 0.2 then the noise makes up less than 20% of the 

profile depth. We rejected spectra with > 0.2. 

3. The placement of the continuum is also an indication of how 
well the CCF was defined. The relative error on the contin- 
uum (^r^) gives a measure of the noise in the continuum 
placement. We rejected all spectra with > 0.1. 

These three criteria identify spectra with poorly-defined 
CCFs. Combined they are the first criteria to be applied to the 
SPC dataset resulting in the rejection of a further 3668 spectra 
from the analysis. Figure [T5b shows the resulting HR diagram 
after these spectra were removed. The hot star mis-classifications 
and overdensity at T^ff ^ 3000 K are diminished. 

Figure[T6h and b show histograms of the SNR per pixel (0.03 
or 0.06 A per pixel) for the full sample of FEROS archived spec- 
tra. The SNR has been calculated within the analysis pipeline 
using sigma clipping on either the predefined continuum re- 
gions in SPB or the pseudo-continuum residual in SPC. The 
highest SNR values do not necessarily mean optimum spectra 
and are identified by the quality flags within the pipeline as re- 
quired. Figure \T6k is the histogram of the spectra that were ac- 
cepted as being of good quality and with a well-defined CCF, 
whereas Figure [T6b is the histogram of the spectra rejected as 
being of poor quality or with poorly-defined CCF. Note that the 
y-axis scale is different between the figures. The histogram in 
Figure [T6h peaks between 50 and 100 SNR. There are signifi- 
cantly fewer spectra with SNR<50 in Figure [T6h compared with 



Figure [T6b. indicating the cleaning thus far has indeed identified 
low SNR spectra. However, as Figure [TOk to d shows, even at 
SNR of 50 the internal error is negligible. 




100 200 300 400 500 600 
SNR per pixel 



200 300 400 500 600 
SNR per pixel 



Fig. 16. Histogram of the measured SNR values of the FEROS 
archived spectra in bins of ASNR - 50: a) Spectra defined as 
good quality and with well-defined CCFs (N=15513); b) Spectra 
rejected as poor quality and/or poorly defined CCFS (N=6038). 
The y-axes are to different scale to better display each sample. 



7.2.2. V,w error 

The next criterion to be applied was the Wrad err or (cry^J which 
had been calculated using the prescription in iTonrv & DavisI 
(.1979.) . Figure [T7] shows histograms of the crv,„^ determined for 
each of the FEROS archived spectra of good quality with a well- 
defined CCF. FigurefTTh is a histogram of Vyad error < 10 kms"' 
in bins of 1 kms"', and Figure [TTb is a histogram of crv,.^^ > 
10 kms ' in bins of 10 kms The majority of the spectra have 
a low cTy,^^, (< 0.5 kms"'), and based on the synthetic spectra 
analysis in Section l4!2l Figure [TT] shows that AW, ad ^ 5 kms"' 
correspond to reasonable variations in the stellar parameters (0). 
However there are a significant number of spectra with cry,^,,, 
greater than 10 kms"' (see Figure [TTb ) which correspond to 
much larger uncertainties in the 6 determination. 




Bin = 50 kms ^ 



10 100 200 300 



Fig. 17. Histogram of measured Yrad uncertainty for each of the 
FEROS archived spectra of good quality with a well-defined 
CCF. a) The distribution of the sample with cry ^^ < 10 kms"'. 
b) The distribution of the sample with crv,,,^ > 10 kms"'. 



Hence we decided to reject all spectra with cry > 
10 kms"'. This resulted in a further 3550 spectra being rejected 
from the SPC dataset. The resulting HR diagram is shown in 
FigurefTSb. Again more of the mis-classifications at hot T^ff, par- 
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Fig. 18. Tgfi vs FWHM of the CCF for a) Giants (log g < 3.5 dex), b) Dwarfs (log g >3.5 dex). The nominal temperature values 
for spectral types AO, FO, GO, KO are shown as cyan lines. A CCF FWHM of 40 kms ' is indicated in red. 



ticular for high log g, have been removed and the overdensity at 
Teff ~ 3000 K is considerably diminished. 

7.2.3. FWHM of V™,, CCF 

Astrophysically, the FWHM of the Yrad CCF can be used to 
calculate the rotational velocity (Vsini) of a star. That partic- 
ular calculation was not carried out here, for the purposes of 
expediency, but the FWHM of the CCF was used to identify 
spectra with broadened spectral features that would not be well 
represented by the synthetic grid. As discussed in Section 15.21 
the synthetic spectra grid was generated with no variations in 
Vsini, assuming all s tars to be slow rotators. The analysis in 
iGazzano et al. I MM found that for a sample of FGK dwarfs 
Vsini < 1 1 kms"' produced good results in the stellar param- 
eter determination by MATISSE. At a resolution of R~26,000 
an upper limit of CCF FWHM = 20 kms ' was set. Due to the 
lower resolution (R~15,000), and hence greater masking of the 
effects of Vsini, of the AMBRE:FEROS analysis we sought to 
relax that criterion. 

The comparison of spectral type with Vsini is used to show 
the increase in the number of fast rotators with hotter spectral 
type for both dwarfs and giants (See Chapter 18, Grav 2005). 
For dwarfs, fast rotators begin to appear at approximately F2 in 
spectral type, for giants at approximately G2. Using the stellar 
parameters derived in SPC and the spectral dataset that remained 
after the application of the crv,^,,, rejection criterion, we replicated 
Spectral Type vs Vsini using Teff and the CCF FWHM as prox- 
ies. The spectra were separated into subsamples of dwarfs and 
giants using the derived values for logg. The comparisons are 
shown in Figure [TSh and b. 

The location of the Spectral Types in temperature are indi- 
cated (cyan lines). There is a clear dense grouping of spectra 
between 4000 and 6000 K with CCF FWHM less than 20 for 
both giants and dwarfs. Both samples have a great deal of scat- 
ter, particularly for T^ff greater than 7000 K. Above the expected 
spectral type of F2 for the giants there is more dense scatter at 
higher CCF FWHMs. For the dwarfs the scatter becomes denser 



around 10000 K. A CCF FWHM of 40 kms ' is indicated in 
both cases (red line). This threshold w as selected as being a fac- 
tor of two greater than the threshold in Gazzano et al. (2010), as 
the resolution of the AMBRE:FEROS synthetic grid is approxi- 
mately a factor of two lower in resolution than that study. Hence 
the astrophysical broadening is masked to at least a factor of 2 
in the convolved FEROS spectra. Above the 40 kms ' threshold 
the distribution of the points is sporadic while below this value 
the distribution is more coherent. Howe ver the majority of th e 
points lie below 20 kms ' in line with iGazzano et alT(l2010l) . 
Ultimately we decided to reject all points with a CCF FWHM 
greater than 40 kms"' . 

The effects of this cleaning on the HR diagram are shown 
in Figure [TStl. A significant number of mis-classifications with 
Teff > 7000 K have been removed as expected from the threshold 
shown in Figure[T8] These rejected stars should be analysed with 
a hotter synthetic spectra grid with variations in Vsini. 

7.2.4. FWHM of the stellar spectrum 

The calculation of the spectral FWHM was necessary in order to 
convolve the observed spectra to the resolution of the synthetic 
grid. However the FWHM values also provide an extra degree 
of information with which to understand the spectral dataset. 
The FWHM measurements were separated into FWHM„,fat, 
FWHMj„„„s, and FWHM,,,^.^,,,,,, classifications (Section [4. 3l l. In 
particular the FWHMj„„„s, classification allowed us to identify 
spectra with broad features, independantly of the FWHM of the 
Wrad CCF For each FWHM,,ro„g spectra, the FWHM of the 
medium strength lines was measured (where possible) as well 
as the FWHM of the identified strong lines. Figure [T9h. b and c 
show the relationships between FWHM,„f.^„„„ and CCF FWHM, 
FWHMmedium and Teff, and FWHM.5,,.0,,^ and Teff respectively for 
the dataset as of the application of the crv„, , criterion. 

In Figure[T9b there is a clear trend of increasing CCF FWHM 
with FWHMmedium- The line at FWHM„,e^,„„, = 0.8 mA intersects 
the trend at the threshold of CCF FWHM - 40 kms"' where the 
trend becomes less dense. This provided the first criterion that 
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= 8 mA, and Tgg = 7000 K are also shown. 



we applied such that any spectra with a measured FWHM„,e^,„„, 
greater than 0.8 mA were rejected. 

Figure [TSfe shows the HR diagram after the removal of spec- 
tra based on this criterion. There is some improvement such that 
the density of the hot mis-classifications is decreased. However 
there are still a significant number of mis-classifications at tem- 
peratures greater than 7000 K. Further investigation showed 
that these remaining mis-classifications were connected to the 
'strong' FWHM classification (FWHM.trong)- 

There was no clear trend between the CCF FWHM and 
the FWHMsirong with which to derive a threshold of rejection. 
Instead a comparison was made to Tgif. Figure [T9b compares the 
Teff with the corresponding FWHM,„f.rf,„,„. The lower boundary 
of possible values of T^ff does increase with FWHM,,,;.^,,,,,, and 
at the derived threshold of FWHM,„erf„„„ =0.8 mA this trend in- 
tersects with a Teff - 7000 K. For greater FWHM,„f^/,„,„ values 
the trend of the lowest Teft disappears into greater scatter, essen- 
tially attributing a temperature value to the FWHM,„(.rf„„„ thresh- 
old. This temperature value was used to infer a FWHM,„erf,„„, 
threshold for the strong spectral lines. 

Figure [T9t shows the measured FWHM.5,rong (where pos- 
sible) with the Teff. There is distinct trend of increasing tem- 
perature with increasing FWHM j,„);ig- The nominal temperature 
threshold at 7000 K intersects with this trend at a FWHM,,,.;,,,^ of 
8 mA. At greater FWHMstrong there is an increase in scatter and 
the trend becomes more dispersed. The greater spread in values 
for spectra with FWHMstrong > 8 mA, and for FWHM,„e^/,„,„ > 
0.8 mA, may indicate a greater uncertainty in the parameterisa- 
tion as the boundary of the giid is approached. 

Based on this we can assume that spectra with a FWHMj„o„g 
not more than 8 mA have stellar parameters that fall within the 
synthetic grid and thus are well-defined. Hence this provided the 
second criterion to apply such that spectra with a FWHMj,ronj 
greater than 8 mA were rejected. Figure [TS? shows that indeed 
the mis-classifications between 7000 K and 8000 K are signifi- 
cantly reduced after the apphcation of this criterion. 



for the rejected sample. It is clearly bimodal with the primary 
peak at 8190+924 K (10068 spectra) and the secondary peak at 
3351+679 K (1903 spectra). The datasets corresponding to the 
key rejection criteria are also indicated: FWHM,„prf,„,„ rejections 
(pink), FWRMstrong rejections (yellow) and CCF FWHM rejec- 
tions (cyan). The spectra rejected on the basis of FWHM,„(.rf,„„, 
and FWHMj,„);ig are concentrated in the two peaks. This argues 
that at least 84% of these spectra with broad features have re- 
turned parameters greater than the 8000 K temperature limit of 
the synthetic spectra grid. 

Figure l20b is the HR diagram of the rejected spectra. The 
majority have Teff > 6000 K, but also with \ogg > 2 dex im- 
plying these are sub-giants and dwarf stars. It is possible that 
these spectra have parameters that lie outside the available syn- 
thetic grid parameter space and MATISSE has compensated by 
varying one or all of the stellar parameters away from the true 
values in order to converge on a solution. Combined with the 
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7.3. Understanding the rejected spectra 

Figure l20l explores the types of stars to which the 11971 re- 
jected spectra correspond. Figures l20t compares Teff with cry^j 



Fig. 20. Comparison of the Tetj of the 11971 rejected spectra 
(black) with a) crv,.„,„ b) logg, c) [M/H] and d) [a/Fe]. The 
datasets for the FWHM,„(.rf,„,„ rejections (pink), FWHMstrong re- 
jections (yellow) and CCF FWHM rejections (cyan) are shown. 
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high number of these spectra that have been measured to have 
broadened spectral features, it is possible that these are in fact 
hot and/or fast rotating stars for which the synthetic grid has not 
been designed. 

In Figure l20h 16% of the rejected sample are located near 
the cool temperature limit of the synthetic grid (3351+679 K). 
Are these truly such cool stars? Figure l20t compares the Teff 
with [M/H]. The cool temperature limit spectra show an inter- 
esting trend of decreasing [M/H] with decreasing temperature. 
For spectra with T^ff < 3000 K the [M/H] extends significantly 
outside the grid range. A possible answer is that these are hot 
or fast rotating stars (based on the FWHM measurements) for 
which MATISSE has converged to a significantly cooler tem- 
perature with significantly decreased metallicity. A similar but 
less dramatic argument can be made for the spectra rejected at 
the hot temperature limit. These lie within the [M/H] range for 
the synthetic grid although many are metal-poor to -3.5 dex. 
Combined with the high gravity values of these spectra the low 
metallicity emphasises the possiblity that these stars have been 
mis-classified due to their parameters lying outside the grid pa- 
rameter space. 

The comparison of [a/Fe] to T^ff in Figure l20h adds weight 
to this, as the cool temperature limit group have very high [a/Fe] 
values. For the hot temperature limit group the [a/Fe] are typi- 
cally within the range of the synthetic grid. 

From this examination of the key variables, the main argu- 
ment that can be drawn as to why ~66.1% of the spectra anal- 
ysed in SPC have been rejected by these very strict criteria is 
that they are most likely hot and/or fast rotating stars that are not 
well-represented by the current synthetic grid. Hence MATISSE 
has converged on solutions which mis-classify the spectra in the 
absence of a grid covering the required parameter range. This 
is not an unreasonable conclusion given that FEROS was de- 
signed and built by researchers within the hot star community 
(Kaufer et al. 1999). A visual inspection was carried out on large 
random samples of the spectra that were rejected based on each 
of the rejection criteria. This inspection supported the appro- 
priate application of the criteria as the rejected spectra corre- 
sponded to noisy, peculiar or strong-lined spectra, and there was 
an over-representation of spectra comprising solely of strong 
balmer lines. However to confirm this conclusion these spec- 
tra can be re-analysed at a later date when a synthetic grid with 
Jeff > 8000 K and a synthetic grid with a range in Vsini values 
will be available. 



7.4. Synthetic grid boundaries 

The above rejection criteria succeeded in cleaning the stellar pa- 
rameter sample without losing the well-defined branches of stel- 
lar evolution that are observed in the HR diagram. Finally, the 
stellar parameters determined near the boundaries of the syn- 
thetic grid are inherently less reliable due to there being less syn- 
thetic spectra present to aid in the definition of the parameters. 
Hence the final rejection criteria that were applied are based on 
the stellar parameters themselves. Taking into account bound- 
ary effects in the learning phase of MATISSE, the limits on the 
accepted parameters are defined as follows: 



3000 < Teff < 7625 

1 < logg < 5 

-3.5 < [M/H] < 1 

-0.4 < [a/Fe] < 0.4 if [M/H] > 0.0 

-0.4 <[a/Fe] < 0.8 if -1.0 < [M/i/] < 

0.0 <[a/Fe] < 0.8 if [M/H] < -I 

Three categories were established using these definitions: 

1. TON: Teff values only within the accepted limits; 

2. TGM: Teff, logg & [M/H] values within accepted limits; 

3. TGMA: T^g, logg, [M/H] & [a/Fe] within accepted Umits. 

The final count of spectra for these three categories were: 
TON = 394 spectra (-292 stars), TGM = 97 spectra (-87 stars), 
TGMA = 6017 spectra (-2780 stars). 

In total 12673 spectra were rejected from the SPC analysis. 
The final set of the stellar parameters accepted for delivery to 
ESO was comprised of 6508 spectra (-3087 stars). For the radial 
velocity analysis 11963 spectra (-4505 stars) were determined 
to have radial velocities with crv„„, < 10 kms"' and these values 
were also accepted for delivery to ESO. The full breakdown of 
the spectral number count is given in Table |8] 



Table 8. Number counts and percentages of spectra accepted 
(Acc.) or rejected (Rej.) during the AMBRE:FEROS analysis. 





No. Spectra 
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% SPC 


No. Stars 


FEROS 


21551 






6285 


Rej. pre-SPC 


2370 


11.0 




1118 


Acc. for SPC 


19181 


89.0 




5500 


Rej. post-SPC 


12673 


58.8 


66.1 


2549 




Accepted for ESO Archive 




^ raci 


11963 


55.5 


62.4 


4505 


TON 


394 


1.8 


2.0 


292 


TGM 


97 


0.5 


0.5 


87 


TGMA 


6017 


27.9 


31.4 


2780 


Total 


6508 


30.2 


33.9 


3087 



7.5. Repeated observations 

The FEROS dataset contains many stars with repeated obser- 
vations. Figure |2T] shows histograms of the number of stars per 
number of observations for the TGMA sample. Figure l2Tt shows 
the histogram for the number of stars with observations less than 
or equal to 10, while Figure l2Tb shows the number of stars with 
observations greater than 10. There is a significant number of 
stars with 2 or more observations, and some stars have been ob- 
served a great many times. The repeated observations were used 
to assess whether the AMBRE analysis derived consistent stellar 
parameters for the same star from different spectra. 

All objects with 2 or more observations were extracted from 
the TGMA sample in order to explore the agreement in parame- 
ters for a single object over a number of observations of similar 
SNR. We also extracted from SIMBAD a list of objects within 
the FEROS dataset (based on a 10" coordinate search) that were 
designated as a binary and/or variable star We found 17 binary 
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Fig. 21. Histogram of the humber of repeat observations in the 
TGMA sample for stars with: a) the number of observations less 
than or equal to 10; b) the number of observations greater than 
10. 
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Stars and 111 variable stars in the final sample which we elim- 
inated from the sample of repeat observations. Hence we re- 
moved those stars with possible astrophysical variations in the 
parameters to be left with variations that are solely due to the 
analysis process. (It must be noted that, as the hst in SIMBAD 
is not exhaustive, it is likely that there are still stars for which 
the parameters vary astrophysically within the repeated obser- 
vations sample.) The repeated observations sample was further 
constrained such that each set of repeats had a mean SNR less 
than 200, as the observations above SNR=200 were sparcely 
sampled. Due to the prior application of the rejection criteria 
there were no spectra present with SNR < 40. Finally, the sam- 
ple of repeats comprised of 584 stars equating to 3438 spectra. 
The maximum number of repeats was 73 for a single star in the 
repeats sample. 

This sample represents the internal error of the AMBRE 
analysis hence we compared it to the internal error (cr,„,) that 
we calculated for each spectra based on the SNR, Wrad and 
Normalisation (Section|5]l. 

Figure|22]replicates the relationships in Figure [TOlbut for the 
repeats sample. For each set of repeats the mean and standard 
deviation were calculated for the SNR and stellar parameters. 
The standard deviation of each of the parameters was used to 
represent the change in the parameters (AO) for each set. The 
mean SNR values were binned (binsize = 25) and for each bin 
the 70''' percentile of the A9 was calculated. These values are 
shown in black. Also included in grey for each parameter are 
the internal errors calculated for the TGMA spectra (cr(0),„,, see 
Section|5]l. 

In Figure l22h the variations in Teff with SNR for the re- 
peats sample is greater than the calculated internal error distri- 
bution and almost traces out the upper limit on the internal error 
However, in general these two independant measurements of the 
internal error are on the same order of magnitude confirming the 
internal consistency within the AMBRE:FEROS analysis. 

As defined in Section |5] the internal errors for log g, [M/H] 
and [a/Fe] are calculated as a piecewise function of SNR and this 
is reflected in the distribution of the internal errors for TGMA as 
shown in Figure l22b to d. The lower limit of the internal er- 
rors for each of these three parameters decreases with increasing 
SNR to the defined constant value at SNR =125. This distribu- 
tion is replicated in the variation of each parameter with SNR 
for the repeats sample. At low SNR there is an offset in log g and 
[M/H] between the two sets of values where the internal errors 
are of greater value. This implies that the internal errors may be 
overestimated at low SNR in light of the repeats sample analy- 



Fig.22. Change in 6 with SNR (in bins of 20) for the repeat 
observations within TGMA (black symbols) compared with the 
o'{ff)int values (Section|5]) for TGMA (grey points), a) The Areff 
that 70% of the repeat sample were less than or equal to at each 
SNR (black squares). The trend with SNR is shown by the black 
lines, b) As for a) but for variations in log g (black stars), c) As 
for a) but for variations in [M/H] (black circles), d) As for a) but 
for variations in [a/Fe] (black diamonds). 



sis. For [(T/Fe] the variations due to the repeats sample trace the 
internal errors very well across the range of SNR values. But 
overall, as for Teff, these two independant measurements of the 
internal error for each parameter are in very good agreement . 

The distribution of the variation of the parameters with SNR 
for the repeats sample in Figure |22] also replicates the analysis 
of the internal error due to SNR that was carried out using the 
synthetic spectra sample as shown in Figure [10] For the lowest 
SNR bin for each parameter the Ad are only 2 or 3 times the 
equivalent value in the synthetic case which is reasonable, and as 
the SNR increases the AO decreases significantly also replicating 
the synthetic case. 

The comparison between the two measures of the internal 
error from the observed spectra and with the synthetic spectra 
analysis shows that the analysis carried out within the pipeline 
is consistent such that closely comparable parameters are found 
for different spectra of the same star, and so for diff'erent stars of 
similar stellar parameters, at a particular SNR value. 

8. The AMBREiFEROS stellar parameters 

The final FEROS stellar parameters are shown in Figure l23lin 
six combinations with which to examine the distribution of the 
parameters. 

Figure l23h shows the Hertzsprung-Russell (HR) diagram of 
the TGM (red) and TGMA (black) samples. The branches of 
stellar evolution are distinct although broad in width with some 
small scatter evident throughout the figure. Figure l2Jb compares 
[M/H] with Teff and the majority of the spectra show near solar 
metallicities (-1.0 < [M/H] < 0.5 dex) consistently across the 
temperatures. 

Figure l23b compares [a/Fe] with T^s and shows an interest- 
ing distribution of low (< 0.0 dex) and high (> 0.2 dex) [a/Fe] at 
low temperatures but an even distribution of solar [a/Fe] across 
all temperatures. 
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Fig. 23. The final FEROS stellar parameters for TGM (red) and TGMA (black) samples for a) The Hertzsmng-Russell diagram; 
b) [M/H] vs Teir; c) [a/Fe] vs T^ff, d) [M/H] vs logg; e) [a/Fe] vs logg; f) [M/H] vs [a/Fe]. 



Figure l23o compares [M/H] with log g and the separation 
between the giants and dwarfs is evident. There is larger scatter 
in [M/H] at low logg (giants). 

Figure l23b compares [a/Fe] with logg and again the scatter 
is larger at low log g (giants). 

Figure l23f compares [a/Fe] with [M/H] for which there is 
a clear trend of enhanced [a/Fe] (~ 0.35 dex) at low [M/H] (< 
-1.0 dex). For [M/H] of solar and above the [a/Fe] decreases to 
solar and depleted values. 

Combining three parameters in one graph provides another 
perspective with which to further explore this dataset. Figurel24l 
shows the HR diagram of the TGM and TGMA samples binned 
by colour in [M/H]. The near solar [M/H] bins contain the ma- 
jority of the sample as expected from Figure |23]across the Main 
Sequence and RGB. The samples within the more metal-poor 
bins lie to the left of the stellar evolution branches which agrees 
with stellar evolutionary tracks. Interestingly the majority of the 
Horizontal Branch (HB) are metal-poor ([M/H] - -1.62 dex). 
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Fig. 24. The HR diagram of the TGM and TGMA FEROS stel- 
lar parameter samples binned by colour in [M/H]. The centre of 
each bin are listed in the legend. 



9. Conclusions & future work 

The stellar parameterisation of the FEROS archive completes the 
first phase of the AMBRE Project. Of the 21551 FEROS object 
spectra that were delivered to OCA, 19181 could be analysed 
for their stellar parameters. The quality flags and error analysis 
resulted in a final total of 6508 spectra (~3087 stars) with stellar 
parameters to be delivered to ESO. Also delivered to ESO are 
the radial velocity values for 1 1963 spectra (~4505 stars) 

Approximately 70% of the spectra were rejected from the 
stellar parameterisation analysis. While 28% were rejected due 
to the quality of the spectra being insufficient for analysis, 42% 
were rejected due to astrophysical reasons. From this analysis 
it seems very likely that hot and/or fast rotating stars are the 
favoured observational object for FEROS. As the current syn- 
thetic grid is defined for stars cooler than 8000 K only then the 
high number of rejections was inevitable. In general the AMBRE 
synthetic spectra grid is currently configured for slow -rotating 
FGK stars. Non-standard stars, such as binaries and chemically 
peculiar stars, within the FEROS sample may still pollute the fi- 
nal accepted dataset as we were unable to identify such stars at 
this stage. 

However, stellar parameter extensions to the synthetic spec- 
tra grid are under development and identification tools are also 
being developed such that a larger range of spectra may be reli- 
ably analysed at a later date. 

In summary the work carried out in the analysis of the 
FEROS archived spectra has resulted in: 

1 . the development of a complex and robust analysis pipeline 
for the determination of stellar parameters of large spectral 
datasets through automated iterative spectral reduction and 
MATISSE analysis; 

2. the establishment of tools with which to exploit spectral at- 
lases and libraries for the comprehensive testing and valida- 
tion of the results from the AMBRE analysis pipeline; 
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3. the determination of stellar parameters (Tes, logg, [M/H|, 
[a/Fe]) for 6508 of 21551 FEROS archived spectra to be 
made available in the ESO Archive, including quahty flags. 

4. the determination of radial velocities (Vrad) for 11963 of 
21551 FEROS archived spectra to be made available in the 
ESO Archive. 

The AMBRE pipeline is currently being used to determine 

stellar parameters for the UVES and HARPS archived spectra, 
the next two phases of analysis in the AMBRE Project. The anal- 
ysis of the FLAMES/GIRAFFE archived spectra is expected to 
commence in early 2012. 
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Table A.l. Description of columns in the table of FEROS stellar parameters delivered to ESO. 



Keyword 


Definition 


Value range 


Null value 


Determination 


DP.ID 


ESO data set identifier 








OBJECT 


Object designation as read in ORIGFILE 








TARGJMAME 


Target designation as read in ORIGFILE 








RAJ2000 


Telescope pointing (right ascension, J2000) 


deg 






DEJ2000 


Telescope pointing (declination, J2000) 


deg 






MJD.OBS 


Start of observation date 


Julian Day 






EXPTIME 


Total integration time 


sec 






SNR 


Signal-to-Noise Ratio as estimated by the pipeline 


0-co 


NaN 




SNRJLAG 


Signal-to-Noise Rado quality flag 


C,R 




C=Crude estimate from SPA* , R=Refined estimate from SPC* 


EXTREME^MISSIONilNEJT^AG 


Detection of extreme emission lines. 


T,F 




T=True: detection therefore no analysis carried out. 










F=False: no detection therefore analysis carried out 


EMISSIONilNEJLAG 


Detection of some emission lines 


T,F 




T=True: some emission lines detected but analysis carried out. 










F=False: no detection therefore analysis carried out 


MEANFWHMilNES 


Mean FWHM of absorption lines around 4500 A 


0-0.33 


NaN 


FWHM measured from spectral features (mA) 


MEANFWHMilNESJLAG 


Flag on the mean FWHM 


T,F 




T=True: FWHM > 0.33 or < 0.1 1. Default FWHM values used 










F=False: FWHM < 0.33, > 0.11 


VRAD 


Stellar radial velocity 


-500 to +500 


NaN 


Units=kms"' 


ERR.VRAD 


Error on the radial velocity 


O-oo 


NaN 


If o-yrad > 10, null value used for afl stellar parameters. Units=kms' 


VRAD.CCFJWHM 


FWHM of the CCF between the spectrum and the binary mask 


0-co 


NaN 


Units=kms"' 


VRADJ^LAG 


Quality flag on the radial velocity analysis 


0,1,2,3,4,5 


-99 


0=Excellent determination. ..5=Poor determination 


TEFF 


Stellar effective temperature (T^e) 


3000-7625 


NaN 


Units=K. Null value used if T^g is outside accepted parameter 




as estimated by the pipeline 






limits or if the spectrum is rejected due to quality flags. 


ERR.INT.TEFF 


Effective temperature internal error 


0-co 


NaN 


Units=K. Square root of quadrature sum of intemal errors 










{0'{Tes)mt.snr, Cr{Teff)int,vrad & c{Tes)mt,norm 


ERR^XT.TEFF 


Effective temperature external error 


120 


NaN 


Units =K. Maximum expected error due to external sources 


LOG.G 


SteUar surface gravity (log g) as estimated by the pipeline 


1-4.9 


NaN 


Units=dex. Null value used if log g is outside accepted parameter 










limits or if the spectrum is rejected due to quality flags. 


ERRJNTXOG.G 


Surface gravity internal error 


O-oo 


NaN 


Units=dex. Square root of quadrature sum of intemal errors 










(a-(logg),„,,j„r, 0"(log g),-„,,vrarf & a"(log g),-„,,„orm 


ERR^XTXOG.G 


Surface gravity external error 


0.2 


NaN 


Units=dex. Maximum expected error due to external sources 


M.H 


Mean metallicity [M/H] as estimated by the pipeline 


0-co 


NaN 


Units=dex. Null value used if [M/H] is outside accepted parameter 










limits or if the spectrum is rejected due to quality flags. 


ERR.INT31.H 


Mean metallicity internal error 


0-co 


NaN 


Units=dex. Square root of quadrature sum of intemal errors 










(o-([M/H]),„,.,„,„ o-([M/H]),„,,„™rf & o-([M/H]),„,,„„™ 


ERR^XTJVIJI 


Mean metallicity external error 


0.1 


NaN 


Units=dex. Maximum expected error due to external sources 


ALPHA 


a-elements over iron enrichment ([a/Fe]) 


-0.4 - 0.4 


NaN 


Units=dex. Null value used if [ff/Fe] is outside accepted parameter 




as estimated by the pipeline 






limits or if the spectrum is rejected due to quality flags. 


ERRJNTjy^PHA 


a-elements over iron enrichment intemal error 


O-oo 


NaN 


Units=dex. Square root of quadrature sum of intemal errors 










(cr([a/Fe]),„,,j„r, o"([a/Fe])i„,,vra<< & o"([a/Fe])i„,,„„™ 


err^xtju.pha 


a-elements over iron enrichment external error 


0.1 


NaN 


Units =dex. Maximum expected error due to extemal sources 


CHI2 


log(x^) of the fit between the observed and the 


O-oo 


NaN 


Goodness of fit between final normalised 




reconstructed synthetic spectrum at the MATISSE parameters 






and final reconstmcted spectra 


CHI2.FLAG 


Quality flag on the fit between the observed and the 


0,1,2 


-99 


0=Good fit...2=Poor fit 



reconstmcted synthetic spectram at the MATISSE parameters 
ORIGFILE ESO filename of the original spectram being analysed 

' = Spectral Processing B 
# = Spectral Processing C 



